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NOTATION 


a  A  constant  for  the  Beattie-bndgeman  equation  of  state* 

(see  Table  I,  Equation  [1.8]) 

A  A  constant  for  the  constant  pressure  heat  capacity  equation  for  the 
ideal  state  (see  Table  1) 

A  A  constant  for  the  Beattie- Bridgeman  equation  of  state  (see  Table  1, 
°  Equation  [1.8]) 


b 


6 


B 


o 


P 

o 


V 

o 


A  constant  for  the  Beattie-Bridgeman  equation  of  state  (see  Table  1, 
Equation  [1.8]) 

A  constant  which  characterizes  the  adiabatic  compression  of  the 
liquid 

A  constant  for  the  constant  pressure  heat  capacity  equation  for  the 
ideal  state  (see  Table  1) 

A  constant  for  the  Beattie-Bridgeman  equation  of  state  (see  Table  1, 
Equation  [1.8]) 

A  constant  for  the  Beattie-Bridgeman  equation  of  state  (see  Table  1, 
Equation  [1.8]) 

Sound  speed  in  the  undisturbed  liquid 

Instantaneous  specific  constant  pressure  heat  capacity  of  the  gas 
inside  the  bubble 

Instantaneous  specific  constant  pressure  heat  capacity  of  an  ideal 
gas  inside  the  bubble 

Instantaneous  specific  constant  volume  heat  capacity  of  the  gas 
inside  the  bubble 

Instantaneous  specific  constant  volume  heat  capacity  of  an  ideal  gas 
inside  the  bubble 

Instantaneous  isentropic  sound  speed  in  the  liquid  at  the  cavity 

wall 


C  A  constant  for  the  constant  pressure  heat  capacity  equation  for  the 

ideal  state  (see  Table  1) 

D  A  constant  for  the  constant  pressure  heat  capacity  equation  for  the 

ideal  state,  (see  Table  1) 

H  Instantaneous  specific  enthalpy  of  the  liquid  at  the  cavity  wall 

n  A  constant  which  characterizes  the  adiabatic  compression  of  the 

liquid 

p  General  liquid  pressure 

p^  Pressure  in  the  undisturbed  liquid;  ambient  pressure 

P  Instantaneous  pressure  of  the  gas  inside  the  sphere 
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Pq  Initial  pressure  of  the  gas  inside  the  sphere 

R  Instantaneous  radius  of  the  imploding  sphere 

R  The  gas  constant 

Rq  Initial  radius  of  the  imploding  sphere 

s  Specific  entropy  of  the  gas  inside  the  bubble 

t  Time 

T  Instantaneous  temperature  of  the  gas  inside  the  bubble 

u  Instantaneous  specific  internal  energy  of  the  gas  inside  the  bubble 

U  Instantaneous  velocity  of  the  bubble  wall 

v  Instantaneous  specific  volume  of  the  gas  inside  the  bubble 

vq  Initial  specific  volume  of  the  gas  inside  the  bubble 

c 

y  Ratio  of  specific  heats,  for  the  gas 

v 

p  Liquid  density 

px  Density  of  the  undisturbed  liquid 
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ABSTRACT 


Two  methods  are  presented  for  calculating  the  in¬ 
stantaneous  pressure,  velocity,  acceleration,  and  radius 
associated  with  the  collapse  of  a  spherical  gas-filled 
cavity  in  an  infinite  compressible  liquid.  One  is  based 
on  the  ideal  gas  law,  the  other  is  based  on  the  Beattie- 
Bridgeman  equation  of  state  for  the  gas  inside  the  cavity. 

In  most  cases  the  latter  assumption  must  be  restricted 
to  relatively  mild  implosions.  The  good  agreement  be¬ 
tween  the  two  methods  serves  to  verify  their  validity. 

Included  are  listings  of  the  two  Fortran  IV  computer 
programs  used  to  obtain  numerical  results  of  the  analyses 
based  on  the  ideal  and  Beattie-Bridgeman  gas  models.  The 
influence  of  several  different  gases,  initial  internal 
gas  pressures,  and  liquids  on  the  collapse  is  studied. 

On  the  basis  of  explanations  of  the  resulting  behavior, 
new  methods  of  producing  similar  behavior  are  discussed. 

ADMINISTRATIVE  INFORMATION 

This  work  wps  funded  under  Deep  Submergence  Systems  Project  Office 
Subproject  S-4607,  Task  11896. 

INTRODUCTION 

Due  to  their  low  density  and  high  compressive  strength  under  hydro¬ 
static  loading,  spherical  glass  shells  have  a  promising  application  in  the 

1  2 

design  of  deep  submergence  vehicles.  ’  However,  due  to  the  nature  of 
fracture  of  glass  shells,  the  underwater  environment  created  by  their 
failure  (implosion)  at  great  depths  is  very  similar  to  that  of  an  underwater 
explosion.  Consequently,  the  effects  of  an  imploding  glass  sphere  on 
neighboring  objects,  especially  other  hollow  glass  spheres,  must  be  given 
careful  consideration  if  hollow  glass  spheres  are  ever  to  be  at  all  suitable 
for  all-depth  vehicles.  Because  the  sequence  of  events  associated  with  the 
failure  of  a  single  glass  sphere  in  a  free  liquid  field  can  be  very  closely 
represented  by  a  gas  bubble  implosion,  the  latter  is  of  primary  interest. 

In  this  paper,  an  extensive  theoretical  investigation  of  the  free-field 
implosion  of  a  spherical  gas-filled  cavity  in  liquid  is  presented.  This 
paper  is  intended  to  supplement  a  previous  paper  on  the  subject.'* 


References  are  listed  on  page  75. 
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EQUATIONS  GOVERNING  BUBBLE  WALL 
MOTION 


THEORY 


4 

Gilmore  has  derived  an  ordinary,  second-order,  nonlinear  differ¬ 
ential  equation  which  relates  the  instantaneous  pressure  of  the  gas  inside 
a  collapsing  (norunigrating)  spherical  cavity  in  an  infinite  compressible 
liquid  to  the  instantaneous  radius  of  the  cavity.  Briefly,  the  equation 
is  obtained  by  employing  the  Kirkwood-Bethe  hypotehsis5  and  basic  fluid 
flow  relations  to  solve  the  spherical  wave  equation.  The  details  of  the 
derivation  can  be  found  in  Gilmore's  report. 

If  R  represents  the  instantaneous  radius  of  the  cavity,  P  the  in¬ 
stantaneous  pressure  of  its  boundary,  and  t,  time,  Gilmore's  equation  is 


R 


1  4R  \  3  /dR\2  /  1  dR  \ 
C  dt  /  +  2  Vdt/  \  '  3C  dt/ 


■'('•if) 


R  dH  /  1  dR\ 

C  dt  V  '  C  dt/ 


[1.1] 


C  is  the  local  instantaneous  isentropic  sound  speed  in  the  liquid  at  the 
cavity  wall. 


C 


c 

00 


(n-l)/2n 


[1.2] 


and  H  is  the  local  instantaneous  specific  enthalpy  of  the  liquid  at  the 
cavity  wall. 


H 


n(P.4B)  r/pB  \(n'1)/n 
(n_l)p« 


[1.3] 


c»«  P«»  and  p»  denoted  the  sound  speed,  pressure,  and  density,  respectively, 
associated  with  the  liquid  when  it  is  in  the  undisturbed  state.  B  and  n 
are  constants  (for  water  B  =»  3,000  atmospheres,  n  —  7)  in  the  formula: 


[1.4] 


which  closely  fits  the  isentropic  compression  curve  for  the  pressure  p  and 
density  p  of  »ny  liquids.  (Except  for  very  large  or  very  small  (cavitation) 
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bubbles,  isentropic  hypothesis  for  the  liquid  is  justified,  because 
the  event  occurs  so  quickly  that  there  is  little  time  for  appreciable 
heat  exchange  to  take  place.) 

For  t  >  0,  the  pressure  P  of  the  liquid  at  the  cavity  wall  will  be 
the  same  as  the  pressure  of  the  gas  inside  the  cavity  provided  the  pressure 
of  the  gas  is  uniform  throughout  the  cavity  and  the  effects  of  surface 
tension  and  viscosity  of  the  liquid  are  negligible. 

The  three  Equations  [1.1]  to  [1.3]  establish  one  differential  re¬ 
lationship  between  R(t)  and  P(t).  But  this  alone  is  not  sufficient  to 
determine  the  behavior  of  the  bubble,  so  another  relationship  is  sought. 

By  assuming  that  the  gas  inside  the  bubble  obeys  some  thermodynamic  equation 
of  state  and  that  the  charge  in  specific  entropy  across  the  bubble  wall  is 
negligible  throughout  the  collapse  (for  the  same  reason  that  an  isentropic 
process  in  the  liquid  was  assumed),  it  is  possible  tc  find  two  independent 

relationships  between  the  pressure,  temperature,  and  specific  volume 
4  3 

(proportional  to  —  nR  )  of  the  gas.  These  two  relationships,  taken  with 
Equations  [1.1]  to  [1.3]  and  appropriate  initial  conditions  from  a  determinate 
system  of  equations  which  can  be  solved  numerically  for  the  instantaneous 
temperature,  pressure,  and  specific  volume  associated  with  the  gas  inside 
the  cavity. 

From  a  computational  standpoint,  the  ideal  gas  law  is  an  advantageous 
choice  for  an  equation  of  state.  When  an  ideal  gas  behaves  isentropical ly , 
the  two  independent  thermodynamic  relations,  i.e.,  the  equation  of  state 
and  the  equation  which  describes  a  zero  entropy  change,  can  be  readily 
combined  to  eliminate  temperature  from  the  calculations,  i.e., 

PvY  «  consta  t  =  P  v  y  [1.5] 

o  o  11 

where  Y  is  the  specifu  heat  ratio,  empirically  determined  for  most  gases, 
and  the  subscript  o  refers  to  some  initial  state,  i'ince  th'-  specific- 
volume  varies  as  the  cube  of  the  radius. 


P  can  be  determined  directly  as  a  function  of  R  by  eliminating  (  be¬ 
tween  Equations  [1.5]  and  [1.6],  yielding 
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/K  v*’ 

••  -  1  \,(mU)  l»*’l 

Iquulion  ( I ,  ’ |  Mini  liqut.  ions  |1.1|  to  |l.3|  constitute  *  »r  t  of  simultaneous 
equations  whose  solution  can  he  obtained  numerically,  Smh  m  Mention  ha* 
in  fact  been  obtained  for  mii  air  bubble  in  water,'  ,f>'  ami  la  eaaHy  ex 
tended  to  other  gases  and  liquids.  Result*  ot  -uich  an  extension  are  pre 
sented  later  in  this  paper. 

A  thermodynamic  equation  of  state  which  is  more  accurate  than  the 

ideal  gas  low  for  gases  at  high  pressures,  but  not  quite  as  simple  to 

apply  in  most  cases  is  the  Beatt ie-Hridgeman  equation  of  state.  In  the 

case  of  nitrogen,  for  example,  the  Beatt ie-Hridgeman  equation  is  accurate 

|| 

in  the  pressure  range  from  one  atmosphere  to  IS, 000  atmospheres,  except 
near  the  critical  point. 

In  this  paper  the  behavior  of  a  spheriial  gas* filled  bubble  in  a 
compressible  liquid  will  be  determined  numerically  by  assuming  that  the 
bubble  wall  obeys  liquations  |l.lj  to  [  1.3 J,  that  the  gas  inside  the  bubble 
obeys  the  Beatt ie-Brldgeman  equation  of  state,  and  that  the  expansion  and 
compression  process  of  the  gas  is  iscntroplc  (reversible  and  adiabatic). 

Letting  v,  P.  and  T  represent  specific  volume,  pressure,  and 
temperature,  respectively,  the  Beatt ie-Brldgeman  equation  is 

Pv2  -  «T  [v  .  »0  (l  -  (•  -  -£j)  -  A„  (l  -X)  II. «l 

where  ft  is  the  gas  constant  (■  0.73032  atm  ft^/mole°R).  The  constants  A^, 

B  ,  a,  b,  and  c,  have  been  empirically  determined  for  a  large  number  of 

o  9  10 

gases.  '  Values  for  some  of  these  gases  are  listed  in  Table  1,  Many 
gases  can  be  uniquely  specified  by  these  five  constants  so  that  the  Beatt ie- 
Bridgeman  equation  represents  a  family  of  equations. 

If,  during  the  collapse,  no  heat  is  exchanged  between  the  gas  and 
the  liquid,  then  each  undergoes  an  isentropic  process  which  det  rmines  a 
mathematical  relationship  between  v,  P,  and  T.  This  relationship  is  inde¬ 
pendent  of  the  equation  of  state.  Thermodynamical ly  speaking,  the  change 
in  the  specific  entropy  of  the  gas  is  zero, 

ds  =  o  (1.9) 
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The  i  pe  v'  1 f  I  o  n)inl  anl  prmnrr  brat  capacity  o  can  bp  found  p* 

plldtly  an  a  dilution  of  I',  v.  and  T  (or  any  (mo  of  thr*p,  kIikp  oi\p  run 

bp  p l la Inal ed  by  menu*  of  the  aquation  of  »tate),  It  happen*  (hut  midi  an 

t>*pre«*  (on  for  c  for  a  h*-  which  obey*  the  Heat  tip  hrldgpntan  pollution 
I*  1 1 

appear*  in  "I'healcai  Proce**  1'rlnciplen ,"  The  detail*  of  the  derivation 
nr*  n I ven  In  Appendix  A. 

.  ,.*>  .  as*  (l .  \  .  "M  ..  ft 

i’  i’  ?  v  j7  svV 

•  t  (*♦».•  £)  fr  •  fv)(1i'T7,5f:;v';v")  1,“1 


whore 


|1.16a] 


•  -  Ao  *  »o  #T 
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•  *A  *  bB  Ht  m  ■ 

o  BB0  ? 
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bcB  R 
.  _ £_ 


11.16b] 


11.16c] 


(1 . 16(1] 


[ 1 . I6e] 


-  B  ft  4 
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2c  B  ft 


y'  -  -  bB  ft  ♦  — f 
o  t3 

2bcB  ft 


[1.16f] 


(1.168) 


[l.lbh] 


and  c  is 
P 

i  .e. , 


the  specific  constant  pre.-.  ■  <re  heat  capacity  of  the  ideal  state. 


c°  «  A  ♦  3t  ♦  Cj2  ♦  Dr 
P 


[1.161] 
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A,  A,  C,  and  ft  a r*»  empl rl cal  ly  determined  constants  *“  which  differ  for 
different  gases,  Values  of  A,  A,  C,  and  ft  are  given  in  Table  I  for  some 
gases.  liquation  (1.15)  with  (1.1ft)  substituted  conatraina  P,  v,  ami  T  to 
lie  on  a  surface  which  character! tea  the  adiabatic  behavior  of  the  gas. 

The  liquations  (1.1ft),  (1.15),  |1.H|,  (1.1)  to  (1.5),  ami  the  re  tat. ion- 
ship  (l.o)  between  v  and  It  taken  simultaneously  form  a  determinate  system 
which  can  be  solved  numerically.  The  equations  are  rewritten  here  for  the 
reader's  convenience. 


R  Jdf  (‘  •  u  $)  *  t  (•  -^3?) 


«  li  ( 

!  ♦  i  R  dH  ( 

1  c  ct /  r  dt  [ 

1  -  c 

iff) 

11-1] 

c  .  c 

“  \  P*,48/ 

n-1 

TrT 

[1.21 

n(P«+B)  r 
H  •  [ 

GW 

n-1 

l  "  - 

(1.3] 

PV2  .  RT  (l  -i)]  (j  -  -Sj)-*„  (l  -i) 


VT  ■ T  (w)„ dp 


c  *  c  + 
p  p 


P  T3  1“  2v2  3v3/ 

*  T(v  *  •.  -(^)5  *  Tfr)  t~  7 


(1.6J 

(1.8) 

[1.15] 

[1.16] 


J. _ 

Y  v+  46 


L J-JL^ 


a  =  RT 


[1.16a] 


Aim* 


H  -  -  ♦  *  ftT  -  4- 

o  o  ,j,2 

[1.16b] 

cB  ft 

Y  -  -  bBrt  ftT  -  — ?— 

oo  t2 

[1.16c] 

bcB  ft 
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[1.16d] 

a'  -  ft 

[1.16m] 
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[1 . 16f ] 

T 

2cB  ft 
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[1.16g] 
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...  0 

6  "  ’  3 
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P 
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THE  INTEGRATION 

In  order  to  see  more  clearly  how  to  solve  the  set 

of  Equations  [1.1] 

to  [1,161],  Imagine  reducing  that  system  by  means  of  substitution  to  three 
equations . 

The  first  equation  can  be  obtained  by  substituting  C  (Equation  [1.2]), 

Ji| 

H  (Equation  [1.3]),  and  ^  (Equation  [1.3]  differentiated  with  respect  to 
time)  into  Equation  [1.1].  The  equation  resulting  from  these  substitutions 
is  a  relationship  (actually  Gilmore's  bubble  wall  equation)  between  R,  R, 

R,  P,  and  P  which  can  be  solved  explicitly  for  R,  i.e., 

R  -  6  (R,  R,  P,  P)  [2.1] 


The  second  equation  is  Equation  [1.8],  the  Beattie-Bridgeman  equation 

3  /„  v3 


of  state,  after  substituting  v 


^  for  v  (the  substitution  v  *  vq 

is  obtained  by  solving  Equation  [1.6]  for  v) .  This  new  relationship  be¬ 
tween  P,  R,  and  T  can  be  solved  explicitly  for  T,  i.e., 

T  -  B(R,P)  [2.2] 


The  third  equation  constrains  the  gas  to  behave  isentropically.  To 
obtain  it,  first  use  Equations  [1.16a]  to  [1.16i]  in  Equation  [1.16]  to 
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gat  c  as  a  function  of  v  and  T.  Noxt  substitute  c  «  c  (v,T)  into 
•  p  P  P 

Equation  [1.15].  Finally,  after  carrying  out  the  indicated  partial  dif¬ 
ferentiation  with  the  aid  of  Equation  [1.8],  replace  v  in  this  new  equation 

by  v0(“g~)  to  8et  *  differential  relationship  between  P,  R,  T,  i.e., 

cp(v(R),T)  dT  -  T  3v^<-T<v^R^)  dP 

where 

dP  -  f  (R,  T,  P)dT  [2.3] 

The  three  Equations  [2.1],  [2.2],  and  [2.3]  represent,  respectively, 
Gilmore's  bubble  wall  equation,  the  Beattie-Bridgeman  equation  of  state, 
and  the  condition  of  zero  entropy  change  in  the  gas  (ds  ■  o) .  The  functions 
G,  B,  and  f  appearing  in  these  equations,  although  very  cumbersome,  are 
known  functions  of  their  respective  variables. 

The  system  of  Equations  [2.1]  to  [2.3]  can  be  further  simplified  by 
eliminating  the  variable  T,  temperature,  between  Equations  [2.2]  and  [2.3], 

T  -  B  (R,P)  [2.2] 


dP  a  f  (R,T,P)dT  [2.3] 

as  follows:  First  replace  T  in  f  (R,T,P)  by  T  =  B  (R,P),  giving 

dP  -  f  (R,B(R,P) ,P)dT 

=  *  (R,P)dT  [2.4] 

Using  Equation  [2.2]  again,  dT  can  be  found  in  terms  of  P,  R,  dP,  and  dR. 
Recalling  the  chain  rule, 


AT  38  A  D  ^  38  4t) 

dr  ‘  Va  dR  +  9P  dP 


[2.5] 


9B  9B 

where  and  are  known  functions  which  uji  be  found  in  terms  of  P  and  R; 
call  these  functions  g  and  h,  respectively. 

Then 

dT  -  g(R,P)dR  ♦  h  (R.P)  dP  [2.6J 

Now  use  Equation  [2.6],  the  value  of  dT,  in  Equation  [2.4]  the  result  being, 
dP  -  g(R,P)  *(R,P)  dR+h  (R,P)  4>(R,P)dP  [2.7] 
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Dividing  Equation  [2.7]  by  dt,  gives  an  expression  for 
^  -  g(R,P)  ♦(R.P)  h(R,P)  ♦(R,P)  ^ 

or  using  the  dot  convention  to  represent  differentiation  with  respect  to 
time, 

P  =  g(R»P)  ♦(R.P)R  +  h(R,P)  f(R,P)P  [2.8] 

Solving  Equation  [2.8]  for  P,  yields, 

[2.9] 

Now  Equations  [2.9]  and  [2.1]  are  a  set  of  simultaneous  ordinary  differential 
equations  in  which  R  and  P  are  the  dependent  variables  and  t  is  the  inde¬ 
pendent  variable. 

Conceptually,  Equations  [2.1]  and  [2.9]  are  much  easier  to  solve  than 
Equations  [1.1]  to  [1.6]  because  the  former  set  of  equations  is  a  more 
compact  representation.  The  actual  solution  of  Equations  [1.1]  to  [1.16i] 
need  not  involve  a  direct  reduction  by  hand  to  the  two  Equations  [2.1]  and 
[2.9],  however.  Instead,  this  reduction  can  be  reserved  for  the  computer, 
but  the  reasoning  behind  such  a  reduction  process  is  necessary  in  order  to 
code  the  solution  of  Equations  [1.1]  to  [1.16i]  for  the  computer. 

The  numerical  integration  procedure,  the  method  of  Hamming,  used  to 
integrate  Equations  [2.1]  and  [2.9],  applies  only  to  systems  of  first  order 
ordinary  differential  equations.  Equations  [2.1]  and  [2.9]  were  therefore 
modified  by  introducing  the  new  variable  U,  where 

R  -  U  [3.1] 

Equations  [2.1]  and  [2.9]  then  become,  respectively, 

U  =  G  (R,U,P,P)  [3.2] 

and 

[3.3] 

By  imposing  appropriate  initial  conditions.  Equations  [3.1]  to  [3.3] 
are  numerically  soluble  by  the  method  of  Hanming.  The  general  application 


;(R,P)  ♦(R,P)U 
,-h(R,P)  ♦(R.P) 


g(R,P)  4>(R,F)R 
l-h(R,P)  ♦(R,P! 
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of  this  method  is  discussed  thoroughly  in  "Mathematical  Methods  for  Digital 

Computers"  °  by  Ralston  and  Wilf  and  is  summarized  briefly  in  a  previous 
3 

paper  by  the  writer. 


INITIAL  CONDITIONS 


When  all  the  thermodynamic  characteristics  of  the  gas  and  liquid  have 
been  determined,  three  initial  conditions,  R(0) ,  U(0),  and  P(0),  are 
required  for  the  solution  of  Equations  [3.1]  to  [3.3].  The  initial  situation 
is  brought  about  by  imagining  that  for  all  time  prior  to  t  *  0  there  exists 
an  infinite  expanse  of  compressible  liquid  uniformly  compressed  to  some 
pressure  p^,  and  that  at  time  t  =  0  there  suddenly  appears  in  this  liquid 

a  nonpulsating  spherical  cavity  of  radius  R(0)  filled  with  some  quantity  of 

gas  under  a  pressure  P  .  Such  an  artificially  conceived  situation  leads  to 
some  physically  untenable  consequences.  For  instance,  if  a  point  in  the 
liquid  is  chosen  such  that  it  lies  on  the  bubble  wall  at  t  =0,  then  the 
pressure  P  at  that  point  is 

P  t  *  0 

o 

P  *  for 

Poc  t  <  0 

By  using  Equations  [1.1]  to  [1.3],  Gilmore  has  shown  that  coincident  with 
the  appearance  of  the  bubble  there  will  be  a  relatively  small  inward  jump 
in  the  velocity  of  the  bubble  wall,  i.e.,  if  an  originally  motionless  gas- 
filled  sphere  is  to  obey  Equations  [1.1]  to  [1.3]  for  all  t  >  0  then  it 
cannot  suddenly  appear  without  having  an  initial  inward  wall  velocity  at 
the  instant  it  does  appear.  The  approximate  value  of  this  velocity  jump, 
R(0+),  obtained  from  Equations  [1.1]  to  [1.3]  is  (see  Gilmore's  report  for 
derivation) 


U(0+) 


R(0+) 


Associated  with  this  jump  is,  of  course,  an  infinite  instantaneous  accel¬ 
eration  of  the  bubble  wall.  In  an  effort  to  avoid  the  initial  infinite 

P  -  p 

acceleration,  one  may  choose  — - -  as  the  initial  condition  for  U(0)  and 
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solve  Equations  [3.1]  to  [3.3]  using  initial  conditions  at  t  *  0  rather 
than  at  t  =  0.  This  is  exactly  the  approach  taken  in  this  report,  i.e., 

U(0)  *  U(0+)  « 

^oo  voo 

THE  EULER I AN  VELOCITY  AND  PRESSURE 
FIELDS  IN  THE  LIQUID 

Provided  the  Euleran  velocity  is  considerably  less  than  the  sound 
speed,  an  approximate  method  can  be  used  to  determine  the  Eulerian  velocity 
and  pressure  at  any  standoff  (given  distance  from  the  center  of  the  bubble) 

4 

in  the  liquid.  The  method  was  developed  by  Gilmore  and  has  been  used  and 

3 

discussed  in  a  previous  paper  by  the  writer. 

RESULTS 

The  equations  appearing  in  the  foregoing  analysis  have  been  coded 
in  Fortran  IV  for  the  IBM  7090  digital  computer  to  determine)  numerically 
the  behavior  of  an  imploding  gas  bubble  in  liquid  both  when  the  gas  obeys 
the  Beattie-Bridgeman  equation  of  state  and  when  the  gas  obeys  the  ideal 
gas  law.  Complete  Fortran  IV  listings  of  computer  programs  based  on  both 
models  can  be  found  in  Appendixes  B  and  C.  Data  input  instructions  are 
included.  Bubble  radius,  velocity,  and  pressure  time  histories  calculated 
from  these  programs  appear  in  Figures  1  to  13.  Implosions  involving  several 
types  of  gases  at  various  ambient  and  initial  internal  pressures  are 
represented. 

A  comparison  between  the  results  obtained  using  the  Beattie- 
Bridgeman  equation  and  those  obtained  using  the  ideal  gas  law  is  made  in 
Table  2.  The  influence  which  the  kind  of  gas  inside  the  bubble  and  its 
initial  pressure  have  upon  the  peak  collapse  pressure  is  summarized  in 
Figures  1  and  2  for  depths  from  100  to  20,000  feet  of  water.  The  kind  of 
liquid  which  implodes  on  the  gas  nlso  influences  the  peak  collapse  pressure 
as  shown  by  Figure  3. 

All  the  results  can  be  extended  to  cases  for  spheres  of  any  radius. 
Suppose  that  at  depth  h  a  solution  exists  for  a  sphere  with  initial  radius 
Rq.  The  radius,  velocity,  acceleration,  and  pressure  are  known  functions 
of  time  at  the  bubble  wall  and  some  standoff  in  the  liquid.  If  the  initial 
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(Text  continued  on  page  27.) 
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PEAK  IN  PSI 


PEAK  IN  PS) 


Figure  2  (Continued) 


Figure  3  -  Peak  Internal  Bubble  Pressure  as  a  Function 
Water  Dipth  Showing  the  Trend  in  n  for  Various  Values 
B  when  P  *  1  atm  and  y  *  1.4 

n  * 


Figure  3a 
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PRESSURE  IN  PS 


Figure  3  (Continued) 


radius  is  multiplied  by  A  =  constant,  then  pressure  and  velocity  will 
remain  the  same  if  radius,  time,  and  standoff  are  multiplied  by  A  and 
acceleration  is  divided  by  A. 


DISCUSSION 

Since  the  Beattie-Bridgeman  equation  accurately  distinguishes  be¬ 
tween  various  types  of  gases,  it  provides  a  means  for  determining  the  in¬ 
fluence  which  the  type  of  gas  originally  inside  the  bubble  has  on  the 
collapse.  In  an  effort  to  discover  the  nature  of  this  influence,  the 
numerical  analysis  employing  the  Beattie-Bridgeman  equation  was  carried 
out  for  implosions  of  gas  bubbles  filled  with  argon,  neon,  helium,  nitrogen, 
ammonia,  methane,  propane,  and  butane.  In  each  case  the  liquid  was  water, 
the  initial  sphere  radius  was  1.5  inches,  and  the  initial  internal  gas 
pressure  was  14.7  psi.  It  can  be  seen  from  Table  1  or  2  that  these  gases 
represent  values  of  Y  ranging  from  1.668  to  1.094.  Results  for  argon, 
nitrogen,  anu  butane  are  shown  in  Figures  4  to  7  for  depths  of  100,  500, 

1000  and  3000  feet  of  water.  Peak  internal  gas  bubble  pressures  are  given 
in  Table  2. 

Unfortunately,  the  extent  to  which  an  analysis  of  this  kind  can  be 
carried  is  seriously  limited.  Although  the  thermodynamic  equations,  the 
Beattie-Bridgeman  equation,  and  the  constant  pressure  heat  capacity  equation 
of  the  ideal  state,  are  representative  of  gases  at  very  high  temperatures 
and  pressures,  these  are  quite  often  not  as  high  as  the  values  reached  in 
the  final  stages  of  gas  bubble  implosions.  The  range  of  applicability  of 
the  thermodynamic  equations  depends  upon  the  constants  given  for  each  par¬ 
ticular  gas  (Table  1).  Nitrogen,  for  which  the  equations  are  applicable 
for  pressures  up  to  15,000  atm,  is  the  exception  rather  than  the  rule.  For 
example,  the  pressures  and  temperatures  developed  inside  a  gas  bubble  during 
the  final  stages  of  collapse  at  a  water  depth  of  500  feet  lie  outside  the 
range  of  applicability  of  the  Beattie-Bridgeman  equation  when  the  bubble 
contains  butane,  propane,  methane,  or  ammonia.  At  a  depth  of  1000  feet 
the  thermodynamic  equations  are  applicable  only  to  the  bubble  containing 
nitrogen.  In  addition  to  not  being  applicable  at  very  high  temperatures  and 
pressure^,  the  Beattie-Bridgeman  equation  does  not  hold  near  the  critical 
point.  In  none  of  the  cases  studied,  however,  was  the  critical  point  reached. 
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Figure  Sb 


Figure  5v 
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Figure  5d 


Figure  S  -  B@«ttie-Bridgeman  Analysis  of  the  Collapse  of  a  1.5 
Inch  Radius  Bubble  Filled  with  Argon/Nitrogen  at  1  Atmosphere 
and  520*  Rankine  and  Immersed  in  Water  at  a  Depth  of  500  Feet 
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Figure  7a 


Figure  7b 


Figure  7c 


Figure  7d 


Figure  7  -  Beattie  Brid^eman  Analysis  of  the  Collapse  of  a  1.5 
Inch  Radius  Bubble  Filled  with  Nitrogen  at  1  Atmosphere  and 
520°  Rankine  and  Immersed  in  Water  at  a  Depth  of  3,000  Feet 


Nevertheless,  even  at  relatively  shallow  depths,  it  is  clear  from 
Figures  4  and  5  and  Table  2  that  the  type  of  gas  upon  which  the  liquid  im¬ 
plodes  significantly  affects  the  pressure  developed  during  the  final  stages 
of  implosion.  Moreover,  at  each  depth,  the  peak  gas  pressure  developed  in¬ 
side  the  bubble  is  a  monotone  decreasing  function  of  the  value  of  Y  for  the 
gas.  In  fact,  this  phenomenon  is  so  well  characterized  by  the  value  of  Y, 
that  the  behavior  of  an  argon  (y  ■  1.668)  gas  bubble  implosion,  described 
by  the  Beattie-Bridgeuian  model,  is  essentially  indistinguishable  from  that 
of  a  neon  (y  ■  1.667)  or  a  helium  (Y  =  1.667)  implosion  under  the  same 
circumstances  even  though  there  are  large  differences  among  the  Beattie- 
Bridgeman  constants  for  these  gases.  Aside  from  the  practical  significance 
this  suggests  that  the  ideal  gas  law  for  adiabatic  behavior  can  be  used  to 
determine  the  influence  of  different  gases  on  the  peak  pressure  of  collapse 
As  can  be  seen  from  Table  2,  the  ideal  gas  law  agrees  quite  well  with  the 
Beattie-Bridgeman  equation  in  describing  the  behavior  cf  gases  inside  im¬ 
ploding  bubbles  in  liquid. 

The  ideal  gas  model  was  used  not  only  to  study  the  same  implosions 
studied  by  means  of  the  Beattie-Bridgeman  model,  but  also  to  extend  the 
results  obtained  with  the  Beattie-Bridgeman  model  at  low  depths  to  greater 
depths.  Excluding  compressibility  charts,  the  ideal  gas  law  is  the  only 
practical  means  of  establishing  an  equation  of  state  for  gases  at  those 
temperatures  and  pressures  developed  during  implosions  at  great  depths. 

The  ideal  gas  model  can  be  used  to  extend  results  to  30,000  feet,  but  it 
should  be  noted  that  these  results  are  purely  hypothetical  for  many  gases. 
Beyond  those  depths  at  which  the  Beattie-Bridgeman  model  (with  constants 
given  in  Table  1)  can  be  applied,  chemical  reactions  such  as  dissociation 
and  ionization  (which  violate  the  condition  of  zero  entropy  change  inside 
the  bubble)  may  be  expected  to  have  significant  effects  on  the  collapse. 

The  information  summarized  in  Figures  1  and  2  is  based  on  ideal  gas  be¬ 
havior.  A  more  detailed  description  of  the  effect  of  initial  internal 
pressure  at  depths  of  1000  and  10,000  feet  is  given  in  Figures  8  and 
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9  for  nitrogen  and  initial  gas  pressures  of  1  and  5  atmospheres.  Likewise, 

for  an  initial  internal  pressure  of  1  atmosphere,  a  more  detailed  description 

of  the  influence  which  the  kind  of  gas  inside  the  bubble  exerts  on  the  peak 

* 

collapse  pressure  is  given  for  argon  and  nitrogen  in  Figures  10  and  11  for 
water  depths  of  1000  and  10,000  feet. 

The  relationship  between  Y  and  the  peak  pressure  developed  in  the 
gas  bubble  during  collapse  can  be  roughly  explained  in  terms  of  the  adia¬ 
batic  compressibility  of  gases  (the  fractional  change  in  volume  of  a  gas 
in  a  reversible  adiabatic  compression).  For  any  gas,  the  adiabatic  com¬ 
pressibility  is  defined  as 


If,  for  simplicity,  ideal  behavior  (and  constant  Y)  is  assumed,  then 
PvY  =  constant 

for  an  adiabatic  process.  From  this 

dv  J_  /—)  J— 

‘  v  dP  =  '  v  \ 3P /  '  PY 

s 

Thus  the  compressibility  of  a  gas  undergoing  a  reversible  adiabatic  process, 
such  as  a  gas  bubble  col  laps  is  inversely  proportional  to  the  value  of  Y 
associated  with  the  gas.  Note  that  it  is  also  inversely  proportional  to 
the  pressure  of  the  gas.  It  follows  that  during  the  implosion  those  gases 
with  relatively  large  values  of  y  are  less  compressible  than  those  with 


★ 

The  liquid  overpressures  plotted  in  Figures  9,  11,  and  13  appear  to 
have  superimposed  upon  them  a  very  sharp  spike  near  their  peaks.  Similar 
sharp  spikes  appear  in  the  corresponding  Hulerian  velocity  plot.  The 
spikes  are  actually  not  spikes,  but  points  at  which  the  pressures  and 

4 

velocities  are  multivalued.  As  Gilmore  has  explained,  these  multivalues 
result  from  the  catching  up  and  overtaking  of  characteristics  with  other 
characteristics  which  originated  earlier  at  the  bubble  wall.  The  different 
speeds  of  propagation  of  the  characteristics  are  due  to  the  changing  sonic 
velocity  of  the  liquid. 
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Figure  10  -  Comparison  of  1000  Foot  Water  Depth  Implosions 
of  Spheres  Filled  to  a  Pressure  of  1  Atmosphere  with  Gases 
Whose  Y-Values  are  1.4  and  1.7 
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Figure  11  - 
of  Spheres 


Comparison  of  10,000  Foot  Water  Depth  Implosions 
Filled  to  a  Pressure  of  1  Atm  with  Gases  Whose 
Y-Values  are  1.4  and  1.7 


relatively  small  values,  at  the  same  pressure.  At  any  given  pressure,  the 
relatively  incompressible  gases  (neon,  argon,  helium)  offer  more  resistance 
to  the  inrushing  liquid  than  do  the  relatively  compressible  gases  (butane, 
propane).  This  added  resistance  is  offered  by  the  relatively  incompressible 
gases  not  only  in  the  final  stages  of  collapse,  but  also  in  the  early  and 
intermediate  stages  when  even  a  very  small  decrease  in  the  velocity  of  the 
inrushing  liquid  can  greatly  reduce  the  momentum,  and  thereby  the  pressure, 
in  the  final  stage.  Because  of  the  nature  of  the  mechanism  described  above, 
the  peak  collapse  pressure  can  be  reduced  by  increasing  the  value  of  y,  or, 
according  to  the  idealized  equation  for  adiabatic  compressibility,  by  in¬ 
creasing  the  initial  internal  pressure  of  the  gas  inside  the  sphere. 

It  has  so  far  been  verified  that  a  small  decrease  in  the  compressi¬ 
bility  of  the  substance  inside  the  bubble  can  bring  about  a  considerable 
decrease  in  the  peak  pressure  of  collapse  in  water.  This  suggests  a  simple 
solution  of  the  sympathetic  implosion  problem  in  which  the  liquid  pressure 
field  generated  by  the  collapse  of  one  glass  sphere  causes  the  failure  and 
subsequent  collapse  of  neighboring  glass  spheres.  A  solution  would  be  to 
fill  the  spheres  with  a  substance  which  is  more  incompressible  than  any 
gas.  At  the  same  time  it  is  desirable  to  obtain  maximum  buoyancy  so  that 
the  density  of  this  substance  should  be  at  least  comparable  to  that  of 
gases.  Unfortunately,  no  such  substance  exists.  All  solids  and  liquids 
are  at  least  several  orders  of  magnitude  denser  than  gases.  Thus,  because 
of  minimum  buoyancy  considerations,  it  is  difficult  to  justify  glass  buoyancy 
spheres  filled  with  anything  other  than  a  gas,  unless  weight  compensation 
is  provided  by  imbedding  the  spheres  in  a  suitable  material  considerably 
lighter  than  water. 

It  is  theoretically  possible  to  decrease  the  peak  pressure  by  a 
method  other  than  that  discussed  above.  This  scheme  utilizes  the  changing 
pressure  and  temperature  inside  the  collapsing  bubble  to  produce  a  chemical 
reaction  involving  gases.  The  gaseous  reaction  products  would  indirectly 
decrease  the  compressibility  of  the  gas  mixture  by  directly  increasing  the 
pressure  in  one,  or  a  combination,  of  the  following  two  ways: 

1.  Heat  may  be  a  product  of  the  reaction.  Since  no  heat  s  exchanged 
between  the  liquid  and  the  gas  mixture  during  collapse  all  of  the  heat 
energy  generated  by  the  reaction  would  go  into  raising  the  temperature  of 
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the  gas  mixture  above  that  to  which  it  would  normally  be  raised  by  com¬ 
pression  alone.  If  the  number  of  moles  of  product  equals  the  number  of 
.moles  of  reactant,  and  if  the  behavior  of  the  mixture  roughly  follows  the 
ideal  gas  law 

Pv  =  nRT 

then  the  increase  in  temperature  must  be  ; ccompanied  by  an  increase  in 
pressure  (and  a  slight  increase  in  volume)  to  a  value  above  that  which  it 
would  assume  if  there  were  no  reaction. 

2.  The  total  number  of  mole,  of  the  reacting  gases  may  be  less  than 
the  total  number  of  moles  of  gases  produced.  The  net  result  would  be  an 
increase  in  the  total  number  of  moles  of  gas  mixture.  Assuming  that  the 
heat  of  reaction  is  very  small,  and  that  the  mixture  does  not  deviate 
significantly  from  ideal  behavior,  i.e.,  if  again 

Pv  =  nRT 

then  the  increase  in  the  total  number  of  moles  of  mixture  must  be  accom¬ 
panied  by  an  increase  in  pressure  (and  a  slight  increase  in  volume)  to  a 
v^lue  above  that  which  it  would  assume  if  there  were  no  reaction. 

During  an  actual  chemical  reaction  in  which  all  reactants  and 
products  are  gases  confined  as  in  the  bubble,  the  liberation  or  absorption 
of  heat  and  the  increase  or  decrease  in  the  total  number  of  moles  of 
mixture  generally  occur  simultaneously  and  tend  to  oppose  each  other  to 
maintain  constant  pressure  as  the  reaction  proceeds.  It  is  unlikely,  how¬ 
ever,  that  exactly  constant  pressure  can  be  maintained.  Inside  an  imploding 
gas  bubble  fhe  net  result  of  a  chemical  reaction  may  serve  either  to  in¬ 
crease  or  decrease  the  pressure  above  or  below  that  whi^h  it  would  normally 
be  in  the  absence  of  a  reaction.  The  reaction  can  then  be  favorable  or 
unfavorable  in  arresting  the  collapse. 

All  possible  types  of  such  chemical  reactions  fall  into  two  cate¬ 
gories:  reactions  which  proceed  immediately  (possibly  explosively'  to  com 
pletion,  and  equilibrium  reactions  in  which  the  extent  to  which  the  reaction 
proceeds  is  usually  determined  by  ne  temperature  and  pressure  of  the 
mixture  of  reactants  and  products.  Jt  would  appear  at  first  that  equilibrium 
reactions  have  a  great  advantage  over  explosive  reactions.  Since  an  ex¬ 
plosion  inside  an  imploding  bubble  would  necessarily  lit  triggered  by  the 
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collapse  of  that  hubhlo,  t  hr  Ihlrihllv  <>t  the  explosion  wmiKI  be  approx i 
mutt' I  y  the  s  <i  me*  regard  less  of  the  depth  ot  implosion.  It  tin*  nagnltudr  of 
that  intensity  were  iet  \  nest  tht*  implosion  of  («  bubble  at  a  watet  depth 
of  It). 000  tort  ,  than  a*  a  100  foot  water  depth  that  maun  Iloilo  of  intensity 
would  possibly  !>«'  mort*  devas  t  at  mg  than  a  simple'  Implosion  without  an  In 
tenia  I  explosion. 

* 

FolUeing  this  reasoning,  numerical  calculat  ions  we're’  madt'  to  de' 

tt’niiini'  rough’ '■  the  behavior  at  a  1000  foot  water  depth  of  collapsing 

bubbles  filled  with  N  ,0  (  and  NO,  in  «iu>mivnl  eipil  I  ibrium  at  a  temperature 

of  r>00°K  tvi  pressures  oi  1/«1,  1,  ami  2  atmospheres.  This  particular 

equilibrium  reaction  was  chosen  in  ause  its  properties  are  well  known. 

About  90  percent  of  the  original  gas  mixture  by  weight  consisted  of  N,0. 

►  ** 

which  is  favored  by  low  temperatures  and  high  pressures.  the  results 
showed  the  reaction  to  be  a  perfect  illustration  of  the  opposing  effects  of 
heat  of  reaction  and  change  of  number  of  moles  of  mixture  discu-sed  pre 
viously.  Since  NO,  is  favored  by  high  temperatures  and  low  pressures  and 
since  the  temperature  and  pressure  of  the  mixture  were  simultaneously  in 
creased  by  the  compression,  the  NO,  was  favored  about  as  much  as  the  NO. 

•*  w 


* 

The  calculations  were  very  similar  to  those  made  here  ‘or  inoit  gases 
except  that  additional  equations  were  necessary  to  determine  the  degree  of 
dissociation  of  N  ,0  j  and  t  v>  cotrect  for  the  entropy  introduced  inside  the 

bubble  by  the  reaction.  The  functional  relationship  between  the  constant 
pressure  equilibrium  constant  for  the  mixture  and  "lie  temperature  of  the 

Hi 

mixture  was  based  on  the  observations  of  Hodenst e in .  ‘  The  kinetics  of  the 
reaction  were  ignored  because  the  equilibrium  establishes  itself  quite 
rapidly.  For  initial  internal  pressures  of  1/J  and  1  atmospln  e,  the  peak 
collapse  pressure  was  well  above  that  which  it  would  normally  be  if  there 
were  no  reaction  and  if  the  Y-value  of  the  gas  were  1.7.  Nevertheless,  at 
an  initial  pressure  of  2  atmospheres,  the  peak  pressure  was  about  \/2  that 
of  an  implosion  involving  an  inert  gas  having  a  i  value  of  1.7  and  initial 
pressure  of  2  atmospheres.  The  validity  of  the  equation  for  equi 1 ibt ium 
constant,  however,  is  questionable  above  pressures  of  2  atmospheres  and 
temperatures  o'  500°F.  Moreover,  N  ,0  ,  which  comprised  most  of  the  mixture 

throughout  the  collapse,  does  not  o’  ty  the  ideal  gas  law  very  well  at  the 
high  temperatures  and  pressures  mentioned  above. 
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(Hil  otl  mint  *•  I  v  ,  hiwi’vi'i  ,  vi*rv  llilli'  N  ,t»  evei'  dl*»v»*  iali*d,  *■'  a  *  i  Mill  fit  Mill 

.  4 

fthltt  in  l  lii1  t  1 1 1*  1 1  «im  wn»  not  <»*hi«*v«»d. 

||ir>  *  a i * u I *1 cd  results  nut  <1 1  iii  which  ill*’  Pt|»i  I  1 1 1*  r  i  nm 

mutton  between  NO,  and  N  ,0  toiil<l  no!  lip  effe*tivelv  utilised  I  o  rtirext 

•  r  1 

the  i  mp  I  o»  i  on ,  may  bp  typical  of  alt  equilibrium  reaction*,  Never!  Iip  Ipm*  , 
the*e  ipiiiatlvp  tPaiiUa  should  not  deter  future  invent i gut  ion  in  till*  Ot 
IP*' l  Ion. 

Hip  remainder  >l  tlip  dii**u*»ton  o  devoted  to  Hip  possible  effect* 
of  immer* iiim  Hip  sphere*  in  liquid*  otlipr  Hunt  water. 

The  tltrpp  parameter*  o^,  H,  an*l  n  appearing  In  Hu  equation 


( 


!'♦  It 

T  Vh 


11.4) 


specify  Hip  liquid  which  implodes  upon  t  h*»  gas,  Theoret  ically ,  tin*  l  sen- 
tropic  sound  speed  In  t ho  undisturbed  liquid  can  Iip  determined  from 
liquat  ion  ( 1 . 4  |  ,  1  .  e.  , 


as  a  function  of  o,„,  B,  and  n  so  that  cv  do^s  not  independent ly  specify  the 
liquid.  The  parameter  o.„  influences  only  tin  period  of  collapse,  not  the 
pressure.  This  can  he  readily  verified  by  substituting  Ap^,  where  A  »  con¬ 
stant,  and  i?  t  for  oti  and  t  respectively  in  liquations  [1.1),  (1.3).  Noting 
that  c  becomes  c  /  fT,  the  substitutions  leave  liquation  (1.1)  unchanged. 

The  influence  which  II  and  n  have  on  the  peak  collapse  pressure  are  summarized 
in  Figure  3  in  which  the  peak  collapse  pressures  are  plotted  as  functions 
of  water  depth  (i.e.,  ambient  pressure  determined  in  every  case  by  multi¬ 
plying  depth  by  the  density  of  water) .  The  influences  of  B  and  n  are  not 
nearly  as  pronounced  as  those  of  y  and  P  (Figures  1  and  2),  the  specific 
heat  ratio  and  initial  internal  pressure  of  the  gas  inside  the  sphere,  re¬ 
spectively.  It  appears,  that  decreases  in  the  values  of  B  and  n  result  in 
decreased  peak  internal  bubble  pressures. 

A  rough  explanation,  similar  to  that  made  for  gases,  can  he  made 
for  the  behavior  summarized  in  Figure  3.  It  is  again  based  on  adiabatic 
compressibility,  this  time  for  the  liquid.  The  adiabatic  compressibility, 
defined  as 
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KZ) 

* 

can  be  1  mmi'il  lately  tuUiiUtiul  from  I  quat  t  on  (1.4)  with  the  help  t  the  »-r> 
lat  |oh*h|p 

v 

i’ 

i'  v 

■I' 

Stme  I  mint  i  >M  |  1 ,  4  |  haw  already  been  evaluated  tor  an  adiabatic  process, 
the  result  t  •» 

l  vlv  I  /  »v  \  _  1 

v  dl'  v  v'p/ 

!* 

Thus  the  compress l bl  lity  of  the  liquid  in  increased  by  a  decrease  in  P  (or 
\\),  b,  it.  in  terms  of  compresslhi  i  ity,  the  behavior  of  the  liquid 
is  opposite  to  tbnt  of  the  jins;  an  i no reuse  in  the  compressibility  of  the 
liquid  is  assoc i nt ed  with  a  decrease  in  peak  collapse  pressure, 

In  order  to  illustrate  tbe  above,  a  comparison  may  be  made  between 
two  liquids  of  equal  density,  one  of  which  is  considerably  more  compressible 
than  the  other.  Shortly  after  the  beginning  of  the  collapse  the  velocity 
of  the  liquid  at  the  bubble  wall  for  a  compressible  liquid  is  about  the 
same  as  that  for  an  incompressible  liquid,  further  away  from  the  wall, 
however,  the  particle  velocity  of  a  compressible  liquid  will  be  less  than 
that  of  an  incompressible  liquid.  Moreover,  at  a  given  time  u  greater 
volume  of  liquid  will  be  in  motion  If  the  liquid  is  incompressible.  In 
fact,  the  liquid  at  a  distance  of  about  cnt  from  the  bubble  center  (t  being 
time  beginning  at  the  instant  of  collapse)  will  be  in  motion;  the  incom¬ 
pressible  liquid  has  a  larger  value  of  ciM  than  does  the  compressible  liquid, 
As  a  result,  the  total  momentum  of  the  inrushing  liquid  will  be  higher  for 
incompressible  liquids  than  for  compressible  liquids  at  the  same  instant  in 
time.  As  time  increases,  the  velocity  of  the  incompressible  liquid  becomes 
increasingly  greater  than  that  of  the  compressible  liquid  at  the  same  dis¬ 
tance  from  the  center  of  the  bubble.  In  addition,  increasingly  more  liquid 
is  set  into  motion.  The  result  is  that  the  total  momentum  of  the  inrushing 
liquid  is  increasingly  greater  for  the  incompressible  liquid  up  to  bubble 
minimum.  The  difference  in  the  final  momenta  of  the  two  liquids  at  the 
bubble  w 3 1 1  influences  the  difference  in  peak  internal  bubble  pressures. 
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keeping  in  mind  the  sympat hefl c  Implosion  problem  mentioned  earlier, 
the  behavior  described  above  suggests  surround i ng  buoyancy  spheres  with  a 
very  compressible  substance.  lino  of  Much  it  substance,  however,  might  lead 
to  a  large  decrease  in  buoyancy  at  great  depths. 

lunation  (1,4)  Is  a  modification  of  the  Tuit  equation  of  state  for  u 
liquid  undergoing  an  isentroplc  process.  The  writer  was  unable  to  find  the 
pa"«meters  B  and  n  tabulated  for  any  liquids  other  than  water.  However,  It 
may  he  possible  to  obtain  good  estimates  of  their  values  for  several  hy¬ 
draulic  liquids  by  fitting  liquation  (1.4)  to  the  compressibility  data  of 

..  . 

Hayward . 

SUMMARY  AND  CONCLUSIONS 

1.  ,A  complete  set  of  equations  has  been  derived  for  the  isentroplc  behavior 
of  u  Beatt ie-Bridgcmun  gas. 

2.  A  Fortran  IV  computer  program  has  been  coded  (see  Appendix  Cl  for  the 
IBM  7090  digital  computer  to  determine  the  behavior  of  a  collapsing  gas 
bubble  in  liquid  when  the  gas  obeys  the  ideul  gas  law. 

3.  A  Fortran  IV  computer  program  has  been  coded  (see  Appendix  B)  for  the 
IBM  7090  digitul  computer  to  determine  the  behavior  of  u  collapsing  gus 
bubble  in  liquid  when  the  gas  obeys  tho  Beatt i e-Bridgeman  equation  of  state. 

4.  Results  of  Items  2  und/or  3  indicate  thut : 

a.  the  ideal  gus  law  provides  a  reasonably  accurate  description  of  the 
gus  inside  a  collapsing  bubble. 

b.  increasing  the  initial  internal  pressure  of  the  gas  inside  the 
bubble  effectively  decreases  the  peak  collapse  pressure. 

c.  increasing  the  value  of  y  of  the  gas  inside  the  bubble  effectively 
decreases  the  peak  collapse  pressure. 

d.  decreasing  the  values  of  B  and  n  in  the  equation  for  isentropic 
compression  of  the  liquid  (liquation  1.4)  decreases  the  peak  collapse  pressure 
somewhat.  Figures  12  and  13  demonstrate  the  extent  to  which  the  peak 
collapse  pressure  can  be  reduced  simply  by  filling  the  spheres  with  argon 

at  10  atmospheres  instead  of  air  at  1  atmosphere.  Detailed  comparisons  are 
made  at  water  depths  of  1000  and  10,000  feet. 
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Figure  12  -  Comparison  Between  1C00  Foot  Water  Depth  Collapses 
of  Spheres  Filled  with  Nitrogen  at  a  Pressure  of  1  Atmosphere 
and  Argon  at  a  Pressure  of  10  Atmospheres 


Figure  13c 


Figure  13d 


Figure  13  -  Comparison  Between  10,000  Foot  Water  Depth  Collapses 
of  Spheres  Filled  with  Nitrogen  at  a  Pressure  of  1  Atmosphere 
and  Argon  at  a  Pressure  of  10  Atmospheres 
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5.  The  results  listed  in  4  verify  th\t  the  peak  pressure  associated  with  a 
gas  bubble  eollupse  in  liquid  can  be  decreased  by  decreasing  the  adiabatic 
compressibility  of  the  gas  inside  the  bubble  and/or  increasing  the  adiabatic 
compressibility  of  the  liquid  in  which  the  bubble  is  immersed. 

t.  Tentative  calculations  indicate  that  chemical  reactions  might  be  utilized 
to  achieve  the  effect  described  in  Item  5. 
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APPENDIX  A 


DETERMINATION  OP  AN  EXPRESSION  FOR  c  «  c  (v,T)  FOR  A 

P  P 

GAS  WHICH  OBEYS  THE  BEATTIE. BRIDGEMAN  EQUATION  OF  STATE 


To  perform  the  derivation  in  an  orderly  manner  it  is  first  necessary 

to  determine  three  fundamental  relations  between  e  and  c  and  P,  v,  and  T. 

P  v  11  14 

These  relations  and  their  derivatives  can  be  found  in  the  references.  * 
Relationship  1:  For  any  gas, 


which  can  be  derived  as  follows.  The  definition  of  c is 


By  differentiating  this  equation  with  respect  to  v  and  holding  T  constant 
there  results 


The  variable  s  is  assumed  to  be  at  least  a  class  II  function  (continuous 
with  continuous  derivatives  up  to  and  including  second  order)  with  respect 
to  the  variables  T  and  v  so  that  the  order  of  differentiation  may  be  inter¬ 
changed. 

Using  the  well-known  Maxwell  relation 


in  the  previous  expression  yields 
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Relationship  2:  For  any  gas, 


as  demonstrated  below.  Assume  that 

s  *  s  (T,v) 


By  the  chain  rule 


ds  *  (ff)v  dT  *  (fr)T  dv 

Dividing  by  dT  and  holding  P  constant  gives 


Multiplying  through  by  T  gives 


This  equation  can  be  rewritten  as 


by  virtue  of  the  definitions 


Finally,  use  of  the  Maxwell  relation 


gives 


48 


Relationship  3:  For  an  ideal  gas 


-  c 


This  can  easily  be  shown  by  employing  Relationship  2  and  the  equation  of 
state  for  an  ideal  gas 

Pv  =  RT 


Substitution  of  the  value 


T 


Pv 

ft 


and  the  derivatives 


into  Relationship  2  yields 


The  superscript  o  indicates  that  the  variables  refer  to  the  ideal  state. 


Now  an  explicit  expression 


can  be  determined  by  sub¬ 


stitution  of  the  Beattie-Bridgeman  equation  of  state  [1.8]  into  Relation¬ 
ship  1,  yielding 


where 


% 


i 
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6cB  R 
0 


6bcB  R 
o 


This  equation  can  be  integrated  from  v  =  °°  (P  =  0)  to  v  =  v  along  a  path  of 
constant  temperature  (note  that  6",  y"»  and  6"  are  all  functions  of  tem¬ 
perature  only) . 


•  v=v 

(ir)T  dv  ‘  ‘  >cv<''-T>lv, 

»  v— 00 


TB"  TY"  T6’ 


2v 


3v3J 


v=v 


V=° 


or 


cv(v,T)  -  [cv(v,T)]v=o 


As  the  pressure  approaches  zero  (volume  approaches  infinity)  the  properties 
of  a  real  gas  become  less  distinguishable  from  those  of  an  ideal  gas  in  the 
same  state;  in  the  limit,  the  corresponding  properties  are  the  same.  If 
c°(T)  is  the  specific  constant  volume  heat  capacity  of  the  ideal  state, 
then 


[r/v.T) 


o 

c 

V 


It  is  understood  that  c°,  by  definition  of  ideal  state,  depends  on  tem¬ 
perature  exclusively.  Then 


Substitution  into  Relationship  3, 

o 

c 

v 


R 
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gives 


This  expression  for  cy  can  be  substituted  into  Relationship  2  after  Re¬ 
lationship  2  has  been  evaluated  for  a  gas  obeying  the  Beattie-Bridgeman 
equation.  The  final  expression  is  the  one  which  appears  in  the  text. 
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AITI  NiUX  « 

COMPHTI'H  I'KOtiHAM  HASH)  ON  111!  IM  A!  Hi  HHIDOIMAN  MODI  l 

Hu*  Fortran  IV  computer  program  HI  10  A  lias  been  coded  to  viol  ermi  no 
numerically  Hit*  behav  i  or  of  a  collapsing  gas  fi  lied  cavity  in  liquid  when 
the  Kas  obeys  the  Heat  t  ie  Hr  l  dgeman  equation  of  state.  The  program  is 
listed  on  the  following  pages.  A  Fortran  IV  or  a  binary  deck  o..»  be  oh 
tamed  from  the  NSRIM!  Applied  Mathematics  laboratory.  A  Fortran  IV  deck  t 

t>e  used  on  any  IHM  /0‘JU  digital  computer  can  he  punched  from  the  listing 

\  8 

provided  the  plotting  routine  gpl^T  is  eliminated  by  following  the  in 
struct  ions  on  four  of  the  comirmt  cards.  Note  that  if  these  instructions 
are  followed,  the  subprogram  CdM,<fT  and  SI’AUN  can  be  excluded. 
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ITERNAl  nfSSURe  IN  THI  SPNCUE  |S,fS.4,ftH  Pftl  A///4BX.2IHANO  THI!  GAS 
4  I  N  S  I  C  E  IS.  I  * , JVS///A0* . PPPAT  AN  INITIAL  TEMPERATURE  Of.fft.O.IGM  0 
<PCRi.es  UNPINS) 
rpr i  r c  to.  i a; 

I  a  r  QRNAT  4  I  HI  /  I  OX,  AliT  IHC.  AX.AP  4  AC  I  os,  PX  ,  BHVtLOC  I  T  ¥  ,  ft  X  ,  I  2MACCELt  WAT  I  ON 
I  .ft  *4  .  Ilf  n»l  SftORC.t.A  ,  I  I  H  f  E  MPE R  Af  ORP  .  l«H  COMPRESS  I  Bit  I  TV. 2* .  I  OH  HIGH  AT 
a  lON/UX  .  I  AM  VIU  I  SI  ClINOS  )  ,  JX.BM  INCHCSI  .  6  X  .  BH I  I  N  /  SEC  )  ,ftX  ,  l«H(  iN/ftKC 
l/Sf l  I .AX  .1 3H4LHG/SOIN) . PX. PHI CEG.R I.AK.VHI SO  IN/LB  I . SX.RH  |  INCHES)/ 
4//) 

C  CALC  WU  UV»|i»0  NPYTUN  TYPE  PROCPOORI  fOR  flNOING  ROOTS  Of  POLYNOMIALS 
VO*  I  4  *  P*  -K.*T  CMOu/PO 
OO  £  1  J-l  ,  ICO 

P  VC 'PC  tv  4  44/14.  P  ~  HG*  T  EMPOA  VCI  •  •  J  ♦  I  BH4D4R  G4BHC  /  T(  MPQ* • 2-OG*  TCMPC4NB 
|  UO)  4  V  C  4  ♦  /  ♦  4  R  C  *  T  (  MP04|?BB*eBB0tRG4HRC4NHB0/TeMP04  4f  “OOA4BOAOI*  V0-RG4 
/PBC*i<P’0*pePii  Tpypu**a 

Of  VI. -4.  •PCX'*  1)4*3/14.  P- 3.04m;  ATtROQ*  VO**  9*1, 0*1  RAAO*RG*OBC  /TEMPO  4* 

1  *-Rv*Tt  NPC*«eem*VO*<  RG*T(Mi>0*Bt-B*PM1B0tRii*BBC*IIBB0/TtMP0**a*BBA*HB 

2  AO  ) 

TVCi'VC-l  VC/OfVO 
!■(  l-9't)  OS, 200. 200 

4  4*9  If  I  AOS  (pi  -TVO)-VU*l  .  Of -4  ) * 0  I  ,  >02*802 
202  VO*TVt 

CO  TO  2P 
200  MR  I  I P(f»«*ll) 

211  P0RRAT4  ltl///5X»  .ISM  ITERATION  flip  VO  OOfft  NOT  CONVERGE) 

GO  TO  9)SJ 
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*  7  COM  IliUI 
401  VITTVC 

5»*. J*  I  I  .01-51  *M0*90RTI  -pn  »  > 

C  CALCULATION  Of  T  Hf  C UP PP I C I  ENT S  UP  THf  TAYLOR  SB  R I  Ml  POH  TMe  INITIAL 
C  V  ALUf  9  CP  IPt  HAPIUS 
OGO* I .0/ VO 

n?«-  (4.  .  7 Jftf.  A  )  *1  N*(  PLAII)  ♦(  l  .0-1  (  PO*B  »/«PL  »6  )  I**  I  tf  N-l  .  01  /KM  I  /  (  I  BA 

i - 1 . j  i*dl  *2 .o*nn » 

CPMJM«I  .  /CGC**  j  *  (  .0*OBC/U'MPQ**  3  I  /OGO* • 2 ♦ BP BO*  (  *.  0*BKC/TS  MP 

io**3-bbbi/cgc-2  .o*»oc*teB*6nMci/ tempo**  3 

cpor  N>R<,*TtMPo/ot.;)**  3*2.  o«(hb!ig*rg«  tempo-bo  ao-bbc*rg/ tempo*  *t  i  /ogc 

I  •  •»*  3,0*  <BBA*PBACl-«UB*H8flC)*RG*TeMP0-BBC*MB8(.l*RG/TENP0**2  I  ,'DC.O*  A.  0* 
2B8C*Bt  B*> B80*RG/TEmP0**2 

CPTO»0.  I6B*I  CPHt  A*CPNC  8*5  •  0*  T  EMPO/B  •  0*CPMCC*  (  9  •  0  ♦  TC  M,-0/  ;> .  0  )  **2*CPH 
ICD*(  U.OMEPPO/O.O  I  •«  J|-RG*e  .0*B8C*RG*DG0/TEMP0**3*I I . 0*BBOC*OGC/2 » 
20 -BbhC*HI*B*nC>0**  8/ 3.0  I  ARC**  2*  TEMPO *0G0*  (  HB60*  I  .0/000-066*6600*000) 

3*11. u«4. 0*EBC*OG(VreMPO**3 )«CPNUM/CPUEN 
PNUM*8.  .  *PC/  <  14 .  7*DG11  )-RG*TPMPU-RG*BBBO*BBC*OGO**2/TP  MP0*«8*t  .  0*RG 

I  *B6P0*(Mu  •  enC*nGO**.T/TBMPO**2-RG*TEMPO*nBBO*BBa*OGO**2*BflAO*flBA»Da 

ao**k* 

POf N*-t MTC/ OGO*HG/OGI> *8.0*RG* fBC/ TEMPO** 3*RG*BBBO-RG*B6nO*BBB*OGO* 
l2.0*flG*T<*BC*flOC*OGO/TEMPO**3-2»'<  *RG*BBBQ*B8C*BB8*0GQ**2/TEMP0**J 
Pt»4*.  I  *HI*PNi  «.'<  RO*POEN  ) 

R3«,>  .0*«2**2*(  i  PL*B  I/I  PO*B:  >•*(  I EN-  I  .0  >/|  t.O*EN  >  |/|  3.0*c  I 
l-l  2.  OT.Itit*  l*PN*l  «>t*M  >*R2/<  3 .0  *C*  I  EN-  I  .0  l*OL*RO  1*1*  I  PI  *B  l/IPO»B  I  I  •• 
31  If N-| .  i  I/I I.O*f N) I -I |PO»B)/l PL*B) l**l | En-1 .0) /( «.0*PN| I  I* 
3(2..7Jt»t  «)*P2/<  3.0*C*OL)*I  (PL4BI/IPUTB  1  I  ••  I  I  EN*  I  .  0  I  / 1  8 . 0*r  N  I  I 
08  ALP  A* I  -«BA*BBAUtBUAO/OGO«PQ/I  14. 7*0G(T*  3) I • I 8 . 0*R3 / 1 OGO* *8*R0 1 ♦ 

I 3.0*P6t-U«R2/IOGU*ROI |/|RG*I  I . O/OGO** 2*BBBO/OGO-BBBO*BB»  I **8 I - I  3. 0* 
8BB4U  *R//  I0GU*R0I*P8/I  I  4 . 7*000**  J  I  ♦’8. 0*PO*Rt/ 1  ‘.4.  T*RO*DGO**3I  l/IRG* 

II  I . J /CGL  *  *2 ♦ OBBU/OGO-BEB  *BBBO  I  I 
02 BL T  *»  t .0*enC*OGO*R2/R0 

ALf ACT  |NB4*BBA0-HBAU/0G0-P0/I  14.7*000**31  l/IRG* I  l . O/DGO* *8 ♦ BBBO/DG 
i  o-njnr*nt-B  1 1 
BE  T AO»- PPC  *OGO 

02  U  MP--C2ALF A/J.O* I -03MFT A/8 .0 “08AL  P A* ALP AO* • 8 /O • 0*1  BE  T  4U4  08BET  A/ 
12.0  ♦  ALP  *0 • •  3  *0.’ Of  TA/2?.0*6PT  AO*  AL  P  40*  *8*0  8ALPA/R. 01/18. 0*S0RT(BET 
2 AC* *2/4. JPCCT ALi»*Lr AO** 3/8 7.0  II  I/I  3.0*1 -BE  I AO/8 .0- ALr 40** 3/ . T. 0*  SO 
.IRT  (MET  AC  «*2/4  .0 ♦ 06  T  AO  •  AL  P  AO*  •  3/. 7.01  )*•(  2.0/3.01  I  ♦  I -D2BE  T A /2 . 0- ALP 
4  AO* *8*0.  4LPA/R,  J-l  B6  T  AO •  02BC  T  A / 2 . 0*  AI.P  AO*  •  3*0BBE  T  A /J  7 . 0*  ALP  AO* *8  *B 
Sf  TAU*0<>ALPA/9.0l/l2.0*S0RT(HFTAO**2/*.0*BETAa*ALPAO**3/»7.0l)  I /| 
63.0*1 -BL I  AC/ 8.0 -ALP  AO** 3/27.0- SORT  I 06 T AO* • 2 / 4 . 0*BE TAO* ALP  A 0 * *  3 /2  7 , 
70  I  >**!»•< /3.0  I  1 

08CNU**I 4,0/1  OGO* *3*R0 1*6. 0/|  CGO**8*RO I *  I BBBOt 8 . 0*BBC / TeMPO** 3  I ♦ 3 . 

1  0  *6U  BL *  I  PRB*2 . G • BBC/ TEMPO**  3  I  / 1 DGO*RQ  11*1 2. 0*R8 l *6. 0*BBC /TEMPO* *4 
2*1 6HB *P OPO -B6BQ/ OGO- I .O/OCO** 8  I • 08Tf MP 

02 COEN* 10,0  *RG*T  EMPU/ 1 OGO*  *3*RO 1*18.0/1 OGO**8*RO 1*1 -8BA0*BBBQ*RG*  T 
IP.  MPO-BBC  «RG/ TEMPO  **2  I  *9.0  '<  0G0*»U  1*1  BB  A*BB  A0-6BB*BBB0*RG*  TPMPO-BBB 

2  0*BHC*RG/TEMP0**2 I  I *1 8 .0*R2 I ♦ I R  G/OGO**  3*2 . 0*RG/DGO**8*  «BBB0*8 . 0*BB 
JC/ TEMPO* *3 l-3.0*MBBO*RC/OGO*l BBH- 8 . 0*BBC /TEMP0**3 1-8.0*608*860*888 
40*RG/ TE  MPO • *3  I *02 TEMP 

02CPT0-O  .368*1 5.0 *CPHC 6/ 9.0*8 .0* TEMPO*CPHCC* < 5. 0/0.0 1 **2«3. 0*CPMCC 
I • TEMPO** J*I5.0/0. 01** 3 l*O2TPMP-|8.0*BBC*RG*OGO/TEMPO**4*l l .0*0800* 
8OG0/2.)-PBBO*BBH*OGO**8/3.0 l*C8TRMP- I  8. 0*«BC*RG*0GO/( TEMPO** 3*80  I ♦ 
31  I .0  *EBHT *C GO/2 .0-86BO*6G6*OGO** 2/3.0 ) *2. 0*R2*3. 0*BOBO*BBC«RG*OGO* 


*  */ /  (  nr  *i  t  kpom  .1 1  •  t  « .  o*pGo*ei!«  -  3 .  o  >*2 .  *2*0004  (  bhhd«  i  .o/ooo-  t| 

ItlfHIHHr,  •0CO)*ll  .0  >,*  .  04BHC  41)00/16 MPH4 4  3)*CPNUM/CPDEN*D2TPMP-3.0*HG  I 

ti •  r  e  »  »r*CGQ/R<i*»  euoo*  :  .o/dgo-hbr4h»«bd*dgo  >*  » i .  0*2 .  o*oc.04bbc/  temp  f 

/u**  j  » *  (.•KUP/cprt  .  i*aG»*2*TeMPn*oG'J/po*  i  i . o/oGQ*BBO*unnn*  S 

ttfK.O)  ♦(  I  .  4  2 .046MC  *0  0(1/ T  CMPO**  3  )  *  CPNUM /CPDE  N*  2 .  0*  R  t- 6  .  0*OHC  *NG*  *2  *0 

toe**; /i  ►pc»*2*«pbbu*i  .G/oGo-PB*'*BBBn*r)an )•  <o*Tepp/Tt hpo** . o*Ra/ ho 

1  »  ♦cpnup/i  °otf Hi, *  *2  *  t  r hpo*oco *i b bbo *  t  .o/ogo-«hr*«bhu*ogo > *  <  i  .0*2.0 

2 •BRC*DGO/TepPO**J I *<  OaCNUM/LPCPN-CPNUXIAOaCOEN.'CPOeN***  } 

DiPnm*..  .0/ (  1 « .  r*o<tu  1  ♦  (  e ,0  *pn/  1  i«.  7*i>gq*ro  )+b ,  o*t*o*eBk)o*BBC*DGr 

14*2/ ( TP  wro* ♦?•«'>» -IB .04MG4PUOU4HHB4BBC4DGO** J/C TPMP04424R0) 46.04RG 

2  4TEMRC*l.  C’HC  4CPH  *l)GU4  *2  /  RQ-6.04n«AO*RHA»OOn4*2/Rn  >4( 2.0*R2  >*  (2.C4RG 
3*8PL'C'*Hi  C4CCa4*.'/T|?MP0**  3-RC.~4,t>*RG*HBRO*BBH*BBC*DGO**  3/ TEMPO* *3 

4  -RG4BBC0  4RB0  4CGII4  42  I  402  T  RHP 

Oi'P|)PN«-|i2CPTO/UGO*C  J.04HG/(UC0*H0)-3.0*CPTO/<0GO4RO»*3.04RG*BBHO4 
I  PBB4  0C.U/  MO-6 .0*HG*0«n0  4«BC4|)GO/(  R04TEMP044  3  )  4  I  2 . 04RG4BBBO * BHO *88 C  • 

2QGQ442/  <  .»C«T  6MPW4  43  |  )4y.04R>'*e.>4RG*BnC  /  TPPPQ4444  (BI1tlO4HBO*DG0442~ 

IPCBn*CGi,- I  .o  I  ♦O.'THHP 

r>3*44  .  I /<  6,0  *RO I *< CPT0  4A . 0 *H 3 *PNUM/POEN ) 

R44  3.04  1  34R2/C*! I PL4fl)/<  P0  4H )  )4*(  (  BN-  I  .  0  )  /  <  2.0»EN )  >-2. 04",*  /  !3.0 

1  4RC  > -I  .’.0  7366  4  >  *  t  N4  (  Pt  4  B  I  /  <  (  C  N-  I  .  0  t  *  OL  •  RO*  2  4 . 0  )  •  <  -  <  t  N-  I  .  0)  /  i£N4|  PL 

2  4  H  )  4  4  (  I  N-l . 0  1/  •  N  »  4<  PO 48  14  4(1  .O/EN ) ) * > . 06Pf 4 6 . 04R 3/C  4  C  I  PL  46  > /(PU4B 

n  I  44  (  t  W*-|  .0  >/<  .',04tN>  >  I  4<  2,0  73664  >/<  C40L44 ,0  I 4PJ4I  (  PL  4B  )  /  I  P04B  )  I  4 
44( ( C  N  4  I .  )  /  I  2 . 0  *  6  N  I | -<2 .07362  4)4P24R2/(C44  24QL )* (PL*U> /( (P04  6*  *1.  ) 

P4 • -44  « I /24.j4C°t (I/HII44  24| PNUN/POEN )*4.0*R24*f444.l 4PNLP/ ( HO* POO N 4 
124.  *4(CPTC*2*.  4R446.04R?  *0?  CP  TU“B . 04CP  TO*R  2*  *2 /RO  >  4  I  32 . 3/ (  I2.0  4R 

20) 4CPI04R24I | .0/PCBN402PnuM-PNUM40>POEN/PDCN44>) 

G  CALCULATION  CP  POUR  INITIAL  VALUES 
*X<  I  )*t  .. 

RNK.UMO 
PNMt.PC 
UKMJ-O.v 
X  X  (.')■  s 

RNH.«Rr4N24xXI2)4424XX(2|4R.I4MXI2)4*24XK(2)4*2*R44XK(2|4*2 
PNM.'.PCt  l]24XK|  2  )  4424KXI  2>4Pl*XXt»)4424XX<  2  )  4  4  24P4  4  X  X  (  2  1442 
UNM2  *2  •  ,4H2*XX(  •'  >43.04H34XX(  2  >44244, 04XXI  2I4R44XXI2  )  *42 
XX<  J  )  *  2  ,  '  43 

<lNM|tH(UR24XX(3>4424XX(  3I4R34XXI  3|4424XX(3>4424H44XX(3)442 
PNMI  ■PQ4P24XXI  3)4424XX<  JI4PJ4XXI  3  >  *•  i! 4 X X  ( 3 ) 4 4 2 *P4 4 X X  I  3 ) 4 4 2 
UNMI  i2..4rT2  4XX())4l.u4R34XX(3>44244.04XX(  3>4R44XX|3)442 
XX ( 4 )» 3.  4$ 

RK»ni'4Pi*  4  X  X  (  4  )  4  4  2  4  X  X  (  4  I  4R34XX  (  4  I4424XXI  4  I44  24R44XXI  4  1442 
PKuPOt P. 4XXI4)4424XX(4|4P34XX<4)4424XX(4I4424P44XX(4I442 
UN*£  .04  7  2  4  X  X  (  4  |  4J.04R34XXI  4  )4  4244.04XXI  4  |4R44XX(  4  )4*2 
KKKKM 

7t»2  VNM.3*VI’4(fiNM3/Rll)44  I 
VLX.'i VP4  IHKM?/R|))44  J 
VNHI «VU4| RKHI/Rl | 443 
VK.VII4  <  RN/BC  14*1 
DVNmJ* J.  4VNP34UNM3/RNM 3 
UVKP2*  3.1/  4VNP24UNP2/RNM2 
UVNN 1*3.  4VNP I 4UN* 1 /RNX I 
OVN»3.‘)4**,*UN/HN 
T  KMJ.X  Te  4PI  V6MJ.  PNH  1  I 
IP( TN43I40I ,401 .402 
401  »RI  Tt!  (6.403  > 

403  FORMAT (  .X///43HITf RAT  ION  PQH  TEMPERATURE  DOCS  NOT  CONVERGE) 
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(.0  To 

402  COM  |  MJi 

!NMi:>*t,  4P|  VN«;.PNM.J  » 

|f|I NM2  ) 401 .*01 >406 
406  CONTI  Ml| 

TNMl  nXTEMPIVNMI  .PK«I  I 
I  F  I TKHI 1401 ,401 .404 

404  CONTINUi 

T  N-KV CMP fVK, PN 1 
IF  I T  N>4 , l ,40 t  .4  •» 

405  CONTINUE 

CA|_  l,  BARPR4  (  ALP  A.  H£  I  A . G AMM A , 0  EL  T A , 0 AL F 4 , OHE  T A , 0G4MMA , DOF L T A  ,  TNMJ  ) 
CPBAR-XCPHCI  VNM3,  ?  NM  J  ,  ALT,  •  HE  1  A  .  GAMMA  ,  OEL  T  A,  DALE  A  . DUE  T A , OG AM M A  ,OOE 
I  L  T  A  ) 

OPNMJ-  XL' PPL’ S  (  VNM.1  ■  DV  NM  3  •  CPBAP  .  AL  P  A  ,  BE  T  A  ,  GAUM  A  ,  DEL  T A , OALF A , OBE T A , OG 
l  AMM  A , OCLLT  A ) 

OUNMJ-XACCEL ( RNMJ.UNM J.PNM3.0PNM  3 ) 

CALL  BAMHR A( ALP  A, HE  T  A, G AMM A , OFL T A , 0 ALF A , DRe T A , D6AMM A , DDEL T A , TNM2 I 
CPBAP- AC  LHC ( VM*_ .  TNM?,  ALFA.RE  T  A,  GAMMA.DEl.TA  ,DALFA  .OHFTA , OG  AM  M A , DDE 
I  L  T  A  I 

0°NM^«  KU°MCS ( VNMJ , OVNM# , CPG AR , ALF A . BF T A , G AMM A , DP L  TA ,0 ALFA , QBE  TA.OG 
I AMM  A , COE  L  T  A  > 

0UNM2»XACCELIRNM2,UNM2,PNM2.0PNM2| 

CALL  BARER A I ALF A.HET A, 0 AMM A, OEL T A , OALF A , OBE T A .DGAMMA . DDEL T A . TNM | ) 
CPBAP-ACPHCI VNMl , TMM I , ALFA.CE T  A , GAMMA , DEL  T  A , O ALF  A . 00  c TA, DGAMMA, OOE 
l  L  T  A  1 

0FNM1-X0PR;  >{ VNMl , OvNMI , CPBAR , Al  F A , OE T A . GAMMA .DEL TA ,D ALF A , OBE T A . OG 
1 AMM A . COt  L  T  4 ) 

DUNM|-X*CCEL( RNN 1  , UNM | , PNM I , OPNM I  1 

CALC  HA  T PR A  I ALFA, BE  T  A, GAMMA .OEL  T  A, OALF A , OBE  T  A , DGAMMA , DOEL  T  A , TN ) 
CPBAR- XL PHCI VN. In, ALF  A , HET  A .GAMMA , DFL  T  A , OALF  A , OBE  T  A , DGAMMA .OOELTA) 
OPN-XCPRLSI VN.OVN, CPBAR, ALFA, eCT A, GAMMA, OEL I  A .OALF A ,OBF TA .DGAMMA .0 
l DELTA) 

OUMX  ACCFL  (RN.UN.PN,  OPN  ) 

IFIKKKK-  I  leot ,6e2.6«2 
6«*2  VV<  1  •  I  l»RNMJ 

Y  V ( J , I  1-RNM2 
YY<  J,  l  ) »  RN  M 1 
YT  (4 . 1  »  =-RN 
YV<  I  .2  I-IKM3 

VV ( Z ,2 ) -UNM? 

Y  Y  I  J  «  2 I’lNMl 

Y  Y  (  4 , 2  1 sUN 
ACC  t  I  I  — u GN  M 3 
ACC<2 1-UUNM2 
ACC( Jl-LUNMI 
ACC  I  4  1-r.UN 
VY( 1 , 31-PNM3 
YY|2, J1-PNM2 
YV13, Jl-PNMl 

Y  Y  (  4 , 3 1 -»PN 
TEMPI l  »=  I N  M  3 
TEMPI  2 J-TNM2 
TEMPI  3  > *  TNM 1 
TEMP 1 4 ) - TN 

AOCOMP ( 1  )»- l . 0 /HNM  3 
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ADCONP1 2 ( MNM2*C»NM2 ) 

AOCONPl  i  )« -UNM| ✓( rnh | *qpn« |  ) 

AOCONP  I  4  PNAOPN  ) 

PNR« HK 
CNRaRN 
PNU«UK 
CNL*UN 
PM>>  PK 
CnP«  Pt, 
k*,*k  »o 
N«4 
L  V*0 

681  RriAH«RN>M44.Q*S/3,0*<2.0*UN-UNNl42.0*UNM21 
oaAP*oNA';*«.o*s/ J,o*<  2.0*dun-cunm  i  *2.0  *  dunk  21 

PBAH»Pnv  34  4 ,04S/J,0*( 2 . 0  40PN-CPNM I ♦ 2, 04OPNM2 1 
KOUNtrO 

589  POOM*PHA('-|  12.0/  121  .0*(PNR-CNR1 
POOV»VG*1*COR/Ru!**3 

P00u*U0A*;-  112.0/121  .0*1  PNU-CNU  ) 

HOOP*  PH  Ail-  112.0/ 121. J*<  PNP-CNP  ) 

tbaw= xre ap ( mcov , modp ) 

IF  (rtiAQ.LC.»j.O)C.O  Tu  401 
OMOD V  *  3  •  4FC1OV4PO0U/MO0R 

CALI  t' ARC  RA<  ALFA.  BEI  A, G APR A  .DEL  I  A  .  CALF  A  .  DSC  TA  «  DGAPMA  .  OD6L  T  A  'PAP  I 

.:::r;r;rr"o"'-o“o”v-w . . . 

OMOOU*XACCKL<  MOUR.MOOU, HOOP. CROUP  ) 

COBP* (9. j^PN-PNHZfJ.OAS*! CMU0P4^.0*0PN-DPMN1 > l/B.O 
COBU» (9.  *UN-UNPJ?4J.O*S*IORl)OU4?.0*OUN-DUN*«l » 1/8.0 
court*  (9.  4BN-PNM243.0*S*(  Mr) DU  ♦  2 . 0  *UN- UNM  111/0.0 
1 F  (  AOS  (PPAH-CORi’l  -  AUS(  CORP  |4  |  .0C- 3  1580.  586 . 58  7 

561.  IF*  AOS  <U  AR-C0RU1-AHS(C0RU1  41  .OC-3  1588, 587.587 

587  R BAH* CORP 49 .0/121  .0  * ( RB AR-CURR 1 
USA  It  =  C0>tU  +  9. 0/1/1  ,041  OH  A  R  — CORO  1 
PBAR*C09n49 .0/121  ,04<  PBAR-CORP 1 
KCUNT*Kf  LNT  4  1 
I F  (4  CJUNT- 300  >589,589.4  50 
450  WRI  IE<»,.  452  1 

452  FC9MArt///SX.4JMIHE  INTEGRATION  ITERATION  OOCS  NOT  CONVERGE! 

00  TO 

5BB  XX(N4|  1  a  XX ( N 1 4S 

Y Y ( N4 1 . l  1  *C  ORR4V • 0/ 121,0*1 RHAR-tORR ) 

y  y<  x4 1 ,  ->  i * c oro 4  1.0/121  .o*<  ubar-corl  i 

VF*V0*1VV(L41, 1 l/ROI*#3 

yy(N41 . » 1-CORP4V. 0/121 .0*1 PBAR-CORP! 

TEPPIAm  1*  X  TEMPI  VF.YYIN41,  J)  ) 

OVF=J.  *vr*YY(N41 .2|/yv(N4|, l | 

|  4*n  ALF  *•«*’*  *•  6AVMA,  del  rA,DALFA,DBETA,OGARRA.OOELYA,TFRP|* 

lfcLTAlCPI  r'  '  Vr,TC*P(  **  1  *  ’  aLFA,HETA,GARRA,OELTA  ,  OALF  A  ,OBE  T  A  .OGAPP  A  .  DO 

1 OPFaXCPRLSl VF.OVF. CPF. ALFA, OETA. GANNA. OELTA.OALFA.OOETA.OGAPNA. ODE 

ACC  1 N4 l  1  * X AC CEL ( Y V 1 N4 1 ,  1  ). VV<  N4 I ,21 ,VY(N4 1 ,3} ,OPF 1 


I 


AOCOPP<  Nt  l)c-YY(Nfl,2)/<YY(N41,l  }*OPF  > 

N«A,»  1 
P NUcUHA < 

CMJ*  CCR'l 
PNR=«GA > 

CN»*CCHX 
PNPiPt AP 
C»,P=CCRP 
RNM1>RN»  ' 

RRH.’iRRK  1 
MAPI  «RK 
RNtrriN, i » 

UNRJ»UNK - 
USM?ilJM| 

UM»1  aUN 
UIS»YY  <  N,2 ) 

PMMtPN  t*i 
PN M2>PM«  I 
PM»1  <Pn 
PN«VY(N. J> 

(jPNo-  3«i'MAP2 
opaa'2»i:pap  i 

DPKM1 *CPN 
OPNaOPF 
DUNM  UCIIKR2 

DUNH 1  *  DO A 
DUA* ACC  I  A  > 

I F ( w-99u 1961  ,961.962 

961  IF!ABS(YY(A,J)-YY(N-1.J))-.U50*YY(N-1. 3) >661.661 ,606 
606  S*S/2.  ,i 

LK*LR* 1 

R*PH*  <0j  ,0*«AAi  ib.OARNMl  ♦ 4 0 , 0 *RNM2 ♦  RNM 3 ) /236 . 0* S / 1 26 . 0*  < - 1 5 . 0* LN4 

190.  :>*UN*1A1*>.6*UNM2> 

R3RH*(lf'.0*RA^l J5.04RNN1 ♦ l 0 h . 0*^NM2 ARNM 3 ) /2ft  6 . 04 S / 1 28 . 0* (-3.04LN 

1  -  54 . 0 *IJN  Y  1*2  7  •  0  *UnM?  1 

04  PH* (bV ,«J*UN4l jft.0*UNPl ♦ AO . 0 *UNM 2 *UNM 3 ) / 256 . 0 ♦ S / 1 2 8 . 0* ( - 1 S . 0 * OUAI 
1  ♦90.0ArUAPl*15,..*OUN**2) 

U3RH»  (  1  2  »  0  RUN  ♦  1  Jft  .  0*IJNM  1  ♦10b.  0*UNM2*UNH  3  )  /236.04S/I28.0*  (  -3. 0*  DUN 
1-54,0  *01*  A  P  1  *27  .  40UNP2  I 

PA PH* (80  .0*PNA1  J5.G*PNMl  ♦  * 0 . 0 *PNM 2 4PNM 3 > / 25 6 . 0* S / I  28 . O*  ' - 1 5 . 0 *OPN 
1499.0*C,'AP14  1S.  j*0PNM2) 

P3PH* {  12,0  «PN41  3  5.0  4PAIMI ♦  1  08 . 0*PNM2 4PAIP 3  ) ✓ 2 56 . 0 ♦ S / 1 28 . 04 1-3.040PN 
1-S4.0*CPAM1+2?.w*DPNM2> 

RNP3-R3P1 
RAR^iRNN l 
RMPl »R4PH 
UNPJ*U3PH 
UNP2*UN« 1 
UNWl  asUAMH 
PAPJ-P  )PK 
PNP2*PN»<  1 
PAIP1  *P4PH 

GO  ro  1j: 

962  NK* 1 

S*XX<NK4| )-XX(NK 1 
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HP* XX (NX *2 1 -XX ( NK ♦ I 1 
HINT  1*1'. 

XINre=YY(NK,l)**3*SM(S*HP)**2*YV<NKE|,l>**3-( 2.Q*S*HP*HP««2  >  *yy  ( 
INK,l)**j-S«*2*YY<NK*2.l|**3)/(2.0-*HP*( SEHP) ) • S* ( S* V Y ( NK *2  » I  )  •  *3- 
•>(S*MPI*r*!M(*l,  I  l**i*HP*YY(NX  ,1  S*HP)|*S«*2 

X  I  NT  3*YY  <  rJK  ♦  1  .  n  *  •  3  *MP  ♦  (  -HP*  *2*  YY  (  NK  ,  t  1  •  •  3*  <  HP  **  2- S*  •  2  }  *  Y  Y  <  NK  ♦  1  ,  1  ) 
l**3*S**'*YY(NK*.'.  t  I  ••  J  )/C  2.0*S*<  S  ♦  HP  )  >* HP  *(  HP*  V  Y  <  NK  .  I  )  *  *i-  1  S+HP)  • 

? YV( NK  *t  , | ) •* JvS* YY( NK*2,  1  I **3 )/<  3.  G*S*  < S»HP  1  1 *HP* *2* X  I NT2 

XPINT lijij 

XP I N  T  2  =  X  INT2/VY(NK*I,  t  )**3 
XPlNT3*XINT3/VY< NK*2. I>4*3 
Z  (  NX  1  *  J  . 

Z(NXe1)=(XPInT1*S*<  I  S*HP  )**2*XP 1NT2- ( 2.0*S*HP+HP**2 )* XPINT I -S**2 • 

1  XPI  NT  3  )/  <2 .0*HP*<  SeHP»  >  *S*(  r.*XP  INT3-<  5*MP  >*XPINT2*HP*  XP  INTI  1/(3. 0* 
2HP •  (  5*ll.»  I  >  *S**2  )*772.6 

7<NXE2)-Z<NKE|)»(XPINT2*HP4<-HP**2*XP|NriE(HP**2-S**2)*XPINT2*S**2 
1  *XP|  NtJ)  /  <2  .(•«&*(  S*HP  )  )  *HP*(  HP*XPlNT  l-{  SEHP  )  •  XP  IN  T  2*  S*  XP  I  NTJ  )  /  1 3 . 0 
2*S*<S*HP  )  )  *HP*  1*772,8 
nk»c 

437  NX  *  NX ♦ I 
S*HP 

HP*  X  X  C NK ♦ 1  )-XX(NK] 

XINT1*XInT2 
XIKT2*Kr  NT  3 

X I  NT  3*  Y Y (  NK  «  1  )  J*HPE(  -HP**;**YY«NK-I.  I  )  *•  3*(  HP*  4  2- S*#2  )  4  Y  Y  I  NX  .  1  )  ** 
1 3*S**2*YY(NK*1 . 1 }♦* 31/( 2.0*5* ( S*HP> )*HP*(MP* VY(NX-| . 1  1 ** 3- <  S4HP) * 
2W(  NX  ,  1  )  **3*S*YY«  NX*  l  ,  I  1  • *3  I /  (  3 . 0 *S* I S+HP  )  )  *  HP*  +2+KINT2 

XPI HT 1 *x  P 1 N  ?  2 
XPINT2SXPINT3 

XPI nT3*X INT  3/YYt  NK+1 . 1  )*«3 

Z ( NK* I  I  =  * { NX  >♦< APINI 2*MP*(-HP**2*XPINT1+(HP**2-5**2}*XPINT2+S**2* 

1  XPI  NT  3  )  /  (2  .0  *S*<  S+HP  1  I  *MP+(hP«XPINT1-<  5 +  HP 1  * XP I  NT 2  + S* XPI  NT3  »  /  »  J,  O* 
?S*(S*HPI |*+P**21*77«:.8 
IF(nK+1  .(  T.999)  it;  Tl)  *37 
DO  9s9  UL*  i  •  N 
959  XXILLI^KXCLL  )*I  .vb] 

PE  AO  (5.731  I0EN.5TF.PCR1  T,  TCP  I  T  ,  VCR  l  T 
73  FORMAT  I  I  I3.4ri5.-j  I 

C  IF  TMI  ,  P«i)CR*K  IS  BEING  REPRODUCED  AND  USED  ON  *  COMPUTER  FACILITY 

C  OTHCR  THAN  THAT  AT  N.S.R.O.C.  ELIMINATE  THE  NEXT  3  CARDS 

GO  TO  (luA.IOS.lij  6),  ICEN 

104  CONTINU- 
106  CONTINUE 

WRI TC (6.  70 1 <  XX(L I . YY(L , 1  1. YYC  L, 2 > . ACC (L 1  .  Y Y C L  .  3  1  .  TE  MP C L  1 . ADCOMP( L I 
1  .2<L  >  ,L= 1  .N  ) 

70  FORK AT (7X.F10.8.F13.S.6E15.6) 

C  IF  T  M I  j  PROGRAM  IS  HEING  REPRODUCE!  AND  USEO  ON  A  COMPUTCR  7ACu.l*» 

C  OTHER  THAN  (HAT  AT  N.S.R.O.C.  ELIMINATE  ALL  CARDS  BEGINNING  WITH  THF 

C  NEXT  UP  T *  .  PUT  NOT  INCLUDING.  TEE  CARO  CONTAINING  STATEMENT  NO.  10V 
IF( I OEx.CC. I )GO  TO  109 

105  CONTINU-. 

OATA  SCALE/OHL  Ii.CAR/ 

XL*  J .0 

CALL  SPACE(XX(M,XR,OX) 
dx«?.o*j» 

Y2MAX* Y V  I  1,2  1 
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y2*i N*rr < 1.2) 

Y  JMA  X*YY  (1,3) 

ZMAx*Z!  1  ) 

7MIn=ZI i ) 

OO  70  (>I.K 

IF  (YY(  l  ,2)  .LT  .Y2MAX  *GO  TO  7tl 
Y2NAX*VY( 1.2) 

70  IF(YY(  I  ,  >)  .GT.Y/MINIGO  TO  79 
V2F(  MYYll  ,2  ) 

79  IF(YY( I . 3) .LT.YJMAX >00  TO  77 
Y3FAX*YY (1,3) 

77  IF(/< I).LT.ZMAX) GO  fO  74 
7X*n«2( I ) 

74  IF u 1 I ) .GT ,ZMIN)CO  ro  729 
ZMlN*ZI I ) 

729  V«V04 < VY ( 1 , 1 )/RO ) 44 J 

IF  (  (  ABL! V-VCRI T  )  .LT.  (  VCR  I T4  > . 03 >  )  .AND.  (  »SP(  YY(  I  . 3! -PCRI  T)  .LT.  !PCRI 
IT  40.05) ) .AND. < AHSI TCR IT- TEMPI  I)  )  .L T . ( TCP  I T4 0 . 05 > ) > MR I TF ( 6 .69 > 

69  FORMATI///5X, 1  2  3H  THE  THERMODYNAMIC  STATE  OF  THE  GAS  MAY  BE  TOC  CL 
IOSE  TC  ThE  CRITICAL  REGION  FOR  I  HE  OF A T T I E-9R I DGEMAN  EQUATION  TO  B 
2F  VALIO  ) 

75  CONTINU.. 

CALL  SPACEIrc.YT.CY) 

READ!  S.76IBT  ITLtil  l  )  ,BT  ITLEI  2)  ,b)ITLE(3).HTITLE(4).BTI  TLE IS)  .B  TITLE 
I (6 ) .HTI TLE ( 7) ,HT ITLhIB ) 

76  FORMAT ( 1  A6) 

REAI1(5.7».)T  ITLEXI  I  )  .  T I  TLEX  (  2  )  •  T  l  TLE  X  (  3  )  .  T  I  TLE  X  (  4  )  ,  T I  **„EX(6>  ,TI  TlEV 
1(1) , T I TL  F V ( 2 ) «  T I TLEV! 3 ) , T 1 TLE V ( 4  >.TI TLEV( 5) 

CALL  GPLCTISCALF . 1 . N . 0 . XL . XR , 0 . - . Y T , OX . D Y ) 

CALL  SPACE! Y2MAX.YT.0YT) 

Y2FINM-Y2MIN 

CALL  SPACE! Y2MI 4.YB.0YB ) 

YB»-YB 

UY»  AH  AX  1  (OYB.DYM 

REA0(S.7G)T|TLEV< 1 ). TITLEVI 2) .T5TL6V! 3 ) . T I TLE V( 4 ) . T I TLE V  I  5 ) 

CALL  GPLCT  (SCALF  .2.N.0.XL.XII,  YB.VT.OX.OY) 

CALL  SPACE! V3MAX.VT.0V) 

RE  AO ! S . 7o ) T I TLEV!  1 ) . T I TL E V ( 2 i . T I TLE V ( 3 ) . T I TLE V ( 4 ) , T I TLE  V ( 5  > 

CALL  GFI  IT<SC*L_, 3.N.O. XL.XP.O.v.YT.OX.OY) 

CALL  SPA  CE I/MAX, YT.OY) 

READ! 5, 76) T  I  TLEV (  I ), T ITLEV( 2) .T (TLEV! 3) . TI TLEV| 4) ,TI TLE  V ( 5 ) 

00  791  N'  L-  l ,N 
791  VVINKL.a )«Z1NXL) 

CALL  GPL  CT  I  3  CALL  .A.'y.O.XL.XR.O.'.YT.OX.OY) 

Y4MIN*!'.  0 
Y4MAX-. 

Y3MIN- 

YSMAXao. 

109  00  99  I  *  I  ,  A 

YX- TV (  I  .  I  )  4YY  C  I, 2)442/2. 04VYI  1.1) /CL 4 ( V Y (  I , 3 ) -PL ) 4 (I . 0- I  V Y ! I  ,  3  >-F L 
I )/l2.040L4C*42 ) ) 42.0736E4 

XX  3*  C  4  4  ./VR4424VYI  I .  I  )4424VYI  I . 2 ) 4 (  l  . 0- Y V C  I . 2  > 44 2 / ( 2 . 04C 4 42 )  ) 

I -C442/Y#  4YV! I , I )•( I .O-VV! | ,2 )/C ) 

YV!  I  ,4 ) 3 YK/I C4STF ) »XX J*YX4  42/(C4  4  34STF4* 2  )  4  (  I.O~VX/(STF  4C  4  *  2 ) ♦ I  XX 3 
I  4  42/ ! 2.  4C4  44 i  ) 4( yx  4  44/<  STF4  4  44C444  )  )  ) 

VY!I  ,5)-CL4!YK/jTF  — YVI  I .  A  )  *  *  2  /  2 . 0  >  /  2 . 0  7  36E  *  *DL  /  <  2 . 04C 4  4 2  >  4  < VX/STF- 
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I YYl I .4  >4 «2/2.0 )«42/2.073bE4 

XX ( I  )=XR (  I  I*(srr-Y»(  I. I  >  )/C*I  1. J-VYI I .2)4 YYl I . I >/<C»STF> )4I .0E3 
IFIYYI  I  .4  )  .LT.Y4VAX)CO  TO  9G 
Y*MX«Y*(  1.4) 

96  1F< YY< I .* > .CT.Y4M I  N  )  GO  TU  97 
Y4XlN*Yr( f .4 ) 

97  IFI  YYC  I  .:»)  ,LT.Y‘5MAX)G0  TO  90 
V5MAX»YV( 1.5) 

9fl  !f(YY(l,S).Ct.Y.i»IMCO  TO  99 
Y5VI N* YY (1.5) 

99  CCkT!NU>- 

C  IF  THIS  PROCRAM  IS  BEING  REPRODUCED  AND  USED  ON  A  COMPUTER  FACILITY 
C  OTHER  T  HAN  TF*T  AT  N.S.R.C.C.  ELIMINATE  ALL  CARDS  BEGINNING  M TH  THE 
C  NEXT  UP  TO.  BUT  NO)  INCLUDING.  TEE  CARO  CONTAINING  STATEMENT  KO.  198 
IF< IOEN.EC.I )G0  TO  108 
CALI.  S*>ACe<XX(N  )  .  XR,  CX  ) 

DX« 1 .0*0* 

CALI-  SPACeiYAMAX.VT.OYT) 

YAM  I N*-Y4M1N 

CALL  SPACElYAMIN.YB.OYB) 

YBs-YB 

Q  Y*  AM  AX  3  ICYT  .0Y(<  ) 

HEAD! o.7 '  ) T I TLE  V  (  1 ) . T ! TLEV <  2 ) . T I TL )  / « 3 ). TITLE VI 4)  , TI TLE VIS) 

CALL  GPi_rT(SCAL«-.4.N.0.XL.X9,  Y8.YT.DX.0Y) 

CALL  SPACE (Y5MAX .YT  «OYT ) 

YSM  t N*-Yt  MIN 

CALL  SPACE (V5MIn.YB.DYB > 

YB«-YC 

DYsAMAX i (CYB.DYT ) 

READ  I  5  «  /<>)  T  |  TLEvI  1  ) . T I TL E V < 2 ) . T I TLE V I  3 > . T I TlF V < 4 )  .TI  TLE  VI  5) 

CALL  GPlCT(3C»Li:.5.n.0.XL.XH,  YH.YT.OX.OY) 

IF< IOFn. ' 0.2 )GC  TO  2 
toe  MRITEIU, l 17)STF 

137  FORMAT (IE|////2^X, 41 MEUL  ER I  An  VELOCITY  ANO  OVERPRESSURE  IN  TME//21 
I  X , 24  H  1ICUI0  AT  A  STANDOFF  OF.F7.2.7H  I NCME S/ / ✓ / I 8X . 4H T I  ME . 9X , I 7HE 
2t-l,E  7  1  An  VELOCITY.  OX.  1  2H0 VERPR C SSURE/ I  3X  .  I  4M(  M  ILL  I  SECONDS  )  .3X.19HII 
INC  Hf S  PER  SiCrNu),4X.  ISH(LeS  MEN  SO  IN)///) 

MR  I  r  C  II.  .139)  <  XXI  1  )  .  YYI  I  .4)  .VYI  1.5  ).  I»l.N) 

139  FCRPATIl  X  .F  -  0  ...  ,flX  .F  12.S.  8X.  E  12 .5) 

2  C  OK  7  | AO 
999  3  T (JP 
tKO 


•IHFTC  BAWDY  NCOECK.StC 

subrout  .  *>t:  oarhrai  alfa. beta,  gamma  ,oel  t  a,  half  a  ,dbe  ta  .ogamma  .odelt  a  , 

I TFMP ) 

COMMUA/ ' F3.PL. c.OL.EN.C. CPHCA.CPHCH, CPHCC.CPHCO.BBA , BBS. BBC. B8 AC 

i ,eehn 

ALFA«  R',4  T  •►P 

OET  A«  -Hit  AC*  ')r><V'l*RG*rLMP-eBC4R5/TEMP4  42 

LAMM  A*B-  «4C0^.C-*t04MMF.O4RC*TEMP-FBC*HBPO4RG/ TFMP44  2 

DELI  A-lruflAfcHCTBt  dCI4RC/TE.<PAA2 

dale A«YG 

PHCT A*H|  M0*RC*2. 0*HHC*rC/TEMP*4  ) 

L'GAMM  Am  -i  EB*f  C0{  4RC4  2.04EBC4  0FBO4RG/TEMP4  4  j 
OOF...  T  A*  -  2 . 0  4#8B*MPC  v  EBB04RC/T  CMP  *  4  3 
RFTllRA 
END 
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t  |  Mh  T  <  Kll  *'HI  I 

rost  i  H.  t 

COKVOMM^uiOUli.m  ,fN.(.  CPMl  A.C**PC«.I  PMf  C  .  C  >*Mf  0  •  BttA  .!**>•(' 

I  .PBi'U 

i  t.pt>»p«v/t  t«  •  7*i»u  t 
IH)  |l*  t  -  »  »  Uni 

t-  T(>*M«tt , ♦  i  v  * nf> n.  »(  i .  i)-nn«/v » i  *ti  bp**  i-  ipava*#/ 1 « •  t mu»»oa  i » « o«»ba/  v 
»  j  >  #n  .*.iu.»«t  t/v*  tv  *m«ntj*4  i  .u-Bttn/v  » i 
«  Di<  t *  mm«  i , o «t«u*i  v i  iO-rofl/v  1 1 •  1 1  •*»* ••#*<■. oaipavaa(?/i*»  * t tin *iit  i 
i  i  .0‘-ri*»/v  i  »  *  r r 
f  If  »*!»«  T  i  *>  P-t  »  f  M.'/OP  1  t  MP 
ir  <  a  -o*t ! .  o ,  to .  »*.* 
to  *Teoi»«t, 
c.0  in  < 

yo  iPtAnttt  OP  1 1»  nf»  i -i  i  pmpai  .it  «a  MfOO*  too.  too 
JOO  TtiO'*»TI  »P 
10  COM  I  Mit 
«>  O0  K  TE  mP<  T  I  •  OS 
b  ULTUHN 
E  NO 


•  im-rc  xct'HCt  Nctrcp.sco 

purl  r  i  or  rcwpc t  v ,  r,  alp  a.mcta,  CAPMA.rtPi  t  a.qai  p a.opp ta .ixiaimra .ttokt  ia 
i  i 

COMMON/  • E/fiC.PL.l'.CL  .GN.  C  .CPMCA  .CPNCB.CPItCC  .CPMCO.Hll*  .BOB  ,BHf  •  Bit  AC 
1  •  BBl  l. 

XCPHL-".  'A  7  7  «t  (‘>MCAtCPHCkAA.OAl  /tl.OtCPPCC  Al  •»  ,  Q  A  T  /  y  ,  Q  I  *  A  #  *  C  PMC  0  *  t  S  . 
t  UAT  /  *>  ,  )  ••  1 1  *«..*  AUBLAMG/  TAA  |A  (  |  .{  ,’V*»M*BIIP|  i  •  0*  V*  *  #  )  -PMRABOPO  /  t  J  •  0  • 

/v«*  .  i  i-  <(.*  i  •(  vt<  »ioo-BBnu*oan/ v)  t(MO/v«*«t  >.  o*HRC*no  / 1  vauaajiaioal 
itA*v«*A»rnpiA*yt*  ituc.ARitAAVA  Ayttiort  I  a  a  v  » /<  al»  A*y**  it  / ,  o»oe  I  A*v**»t 
tl.QtUAtMtVM.OtCUIAI 

weruiiR 

LND 


SI0PTC  KCnrat  NCOfcC* ,  S  C  U 

EORt  T  |  ON  RC“He  St  y  .OV.  CP.  At  Pa.  «f  TA.CAOMA.UPL  TA.OALPA.  OflP  ?A  .O&AMMA  ,D 
lOEtl  *  1 

C  00  MON / A  1  /  no. PL .P.DL .PN.C.CPHCA . CPECR • CPHCC . CPMCO ,«U A , UOH .HOC .BRAC 
t  .BHl'O 
T  *  AL  r  A/  I  j 

Pa  1  A  .  7/y  AAS  A  C  «<.»t  *t  V  tBBHOAl  I  *  O-HHH/V  1  »  *1  I  .  Q-BtlC  /  1  y*  f  ♦♦  J  1  »  -BHrO*  t  I  • 

I  0-boa/ V )  ) 

Tpy»>T«rt(v*pppnt(t.i-BeB/v)>*ti  .0  t».o*unc/t  v*t**3)  »/i/.o*p»v/i»a.7 
i  t  ro  1 1  on  a  t  *Etj.A/ 1  n  j*v**2  >-r-eMnn*HPO*i/v**?-nnc*Pmto/|  v*i  )  a a>* 4.0 a pm 
2C tarn  *rtHO/ ivaaiataa*) i 

Coer  1  *  1 1.  ALP  A  tops  T  A/y  trOAPMA/y  ♦•/  tOOPL  TA/y**  .1  »•  TPVOT 

cceF/^c.'tv 

RPS* I  A . / *CF*OVA| alp A/ytJ .OAPF  TA/VAA2*  3 . 0*0 AMM A / V • • Jt  A. 0»OCL I A/yA *A 
I  ) 

X0PPP.S*.1«-5/l  CCP7  »  -«  UEPf  > 

UPTOWN 
f  NC 
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•  |  III  It  Oil* 

PMNttlCX  KAti  PL  IM.W.P.  OH  I 

MIAMI  l*/  *l  ’  R  A  i  PL  t  M .  PI  tit)  t .CPMt  A.CPHCM.CPMC  I  it  PHI  It.  MM  A  iSMR iMMC  *itt  AC 
I  ,  PRlMI 

|l|Mi»l,l(intMI  «IIP/HL*I  I  PI  »M  l/IP.P  I  »  .«'/**«» 

F  MM-  (  f  i  MM<  |*IN*(PL  *W  |/(<  »N*>  I  •  ft  |*tiL  1*1  I  t  H»L  ♦nn  •••<!!«'•  I  *U  I 

»'PM-».  > 

t  f  .(  *  |  )  •  *  p  (/  |  Pi  •KI|tt|||i,'-liOUIflOttNl) 

■  Aft,  PL  -  ■  >  ,«»t!« •f)«n  ,0-U/  .  I.  >»C«r  »  I  /!*••«»  .ll-Al/tl  )  |  |  .0»V/Cd  *  I 

|  M  •  <  I  *0  '/f  C  I  I  ♦  OF  M I H/ t  t 

MF  ItlXF 
»HO 


•  I  ML  ON  6PLIII 

•  I  ML  OX  AMI  CIV 
1 1  ML  OX  APXMV 

•  MUCH  RIILPC 

•  IIILOX  1*111  1C 
I 1 IH.  CM  Vt'HARV 

•  I  ML  OX  HU  IV 
SI  PL  OR  Fk  l)T  y 


GfUT  binary 
d#ek 


0*/l 0/»A 
CA/l M/AM 
OA/IM/AA 


•  I  WP  I L  SPACE!.  ACCICK.SCC 

SUSXIHJ  11*1  IM(H  CNLPI  ,  KR,C*  I 
IFIIMPI  .Ll  .1  .Ot-JIbO  IU  «• 

*IO**LCC  lOI*M'PT*IOt  9.01-J.*' 
ip  <*  ic. li  .o.oiuh  in  ii« 
ir»*  ic •  i  .o 

ftl  TO  III 

it*  i d* *  I r 

IP  t IP.CF . I Ibl  III  JAM 
M|  I  RR.t>pPT*lti.Q**<  l-IPi 
KM* I *  R»  I 

KH.KM/IU  ,0«*t  l -  10  | 

IMI  IHiLt  ill  IOC  It)  »6I 
OF. | KR/jn 

o  h  >o  */ 1 v ,a«*( i-ioi 
tie  io  mi 

461  CK.| .0/1  <,Q»*I I -10) 
lit  1  0  I  11 

46t>  I  *R«  VA(PI/  |v)  ,0«*  I  1C- I  I 
DM* | INI  1 

KR.«R«I  .0 • *  «  IC-  I  I 

ip  ( i  *r. li  .linen  it*  no 

OK. | XR/»  » 

UK.UK.lv  .0**1  IC-  1  I 

ec  h  iii 

I  16  OK. I  .0*10. Q**<  m-l  I 
bC  I U  III 
AM  »R I  I r  t  6 i  i  ll 

61  PCRMAII  IM///I0*.«6H  IHt  iMJPLI  I  IS  TOO  SPALL  POM  THIS  SOBROL  T I NB 
I  I 

III  RPT'JR*. 
t  KD 
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l  or  each  Mil* le  collapse,  informal  Ion  mom  he  read  In  on  data  card* 


in  the  foil  owing  way: 

laru 

Cola 

X 

1  -  to 

collapse  depth  in  feet  of  water, 

F10.4 

Cola 

11-20 

initial  sphere  radius  In  inches, 

1*10. 4 

Cola 

21-50 

initial  internal  gaa  pressure  in  pal. 

P10.4 

Cola 

51-40 

Initial  temperature  in  degrees  Ranklna, 

FI0.4 

Cola 

41.50 

% 

tha  Beattie-Brldgeman  constant  A  In  USS-Xl-t 

0  mole** 

F10.4 

Cola 

51-00 

ft^ 

tha  Baattla-Brldgaman  conatant  Bq  in 

F10.4 

Cola 

01-70 

ft* 

tha  Baattla-Brldgaman  constant  a  in 

F10.4 

Cola 

71-10 

ft* 

tha  Baattla-Brldgaman  conatant  b  in 

F10.4 

Card  2 

Cola  1.15 

ft*  R* 

tha  Baattla-Brldgaman  conatant  c  in  0K  , 

HIS. 4 

Cola 

10-25 

A,  tha  firat  conatant  in  tha  idaal  constant 
praasure  haat  capacity  aquation, 

F10.4 

Cola 

20-40 

1,  tha  aecond  conatant  in  tha  idaal  conatant 
praaaura  haat  capacity  aquation, 

E15.4 

Cola 

41-55 

C,  tha  third  conatant  in  tha  idaal  conatant 
pressure  heat  capacity  equation, 

K1S.4 

Cola 

56-70 

0,  the  fourth  constant  in  the  ideal  conatant 
pressure  heat  capacity  equation, 

E15.4 

Cola 

71-80 

the  name  of  the  gas  inside  the  sphere. 

Card 

JL 

Cola 

1-  9 

blank. 

Col 

10 

1,  2,  or  3  if  plotting  routine  is  incorporated  in 
program;  otherwise,  blank. 

1  for  printed  output  only. 

2  for  plots  only, 

3  for  printed  and  plotted  output. 
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Cola 

11-2 S 

tha  standoff  to  inches, 

IIS, 5 

Col. 

26-40 

the  critical  pressure  of  the  gas  in  pal, 

H15.S 

Cola 

41-55 

the  critical  teaiperature  of  the  gas  in 
degrees  Rankine, 

m.5 

Cola 

56-/0 

ft5 

the  critical  volume  of  the  gaa  in 

HIS. 5 

Cards  4  through  10  contain  graph  lahala.  If  no  plots  ara  desired 

(i.e.,  if  the  number  appearing  in  Col  10  of  Card  2  it  1)  or  if 
the  program  it  not  bring  used  at  N8R0C,  than  these  cardt  mutt  not 
ha  included  in  the  data. 

Card  4 


Cola 

1-46 

Main  graph  title  for  all  graphs. 

Card 

_5 

Cola 

1-30 

Horiaontal  graph  label  for  all  graphs 
(time  in  milliseconds) , 

Cola 

31-60 

Vertical  label  for  radius-time  curve 
(radius  in  inches) . 

Card 

_6_ 

Cola 

1-30 

Vertical  label  for  wall  velocity-time  curve 
(wall  velocity  in  inches  per  second). 

Card 

2_ 

Cola 

1-30 

Vertical  label  for  bubble  wall  pressure-time  curve 
(wall  pressure  in  psi). 

Card 

JL 

Cols 

1-30 

Vertical  label  for  migration-time  curve 
(migration  in  inches). 

Card 

JL 

Cols 

1-30 

Vertical  label  for  Eulerian  velocity-time  curve  at  the 
standoff  given  on  Card  3  in  Cols  11-25 
(Eulerian  velocity  in  in/tec). 

Card 

w 

Cols 

1-30 

Vertical  label  for  overpressure-time  curve  at  the 
standoff  given  on  Card  3  in  Cols  11-25 

(overpreasure  in  psi). 

65 


Two  MenK  certl*  in  #ucce##lon  atop  the  computer,  linch  c«»e  require#  no 
more  then  4  minute#  running  time. 


APPENDIX  C 

COMPUTER  PROGRAM  BASED  ON  THE  IDEAL  GAS  LAW 


The  computer  program  RU02  has  been  coded  in  Fortran  IV  to  determine 
numerically  tha  bahavior  of  a  collapsing  gas  fillad  cavity  in  liquid  whan 
tha  gas  obays  tha  idaal  gas  law.  Tha  program  is  listad  on  tha  following 
pagas.  A  Fortran  IV  or  a  binary  dack  is  obtainabla  in  tha  sane  way  as  a 
dack  for  RU03;  aither  through  tha  NSRDC  Applied  Mathematics  Laboratory,  or 
by  means  of  tha  listing  and  Instruction  coament  cards  therein.  If  tha 
program,  as  listad,  is  run  on  tha  NSRDC  computer  facility,  than  no  option 
between  printed  and  plotted  output  is  available.  Output  is  always  both 
printed  and  plotted  because  tha  plotted  output  may  be  in  error.  Tha 
numbers  labelling  tha  vertical  axis  can  have  no  more  than  6  digits,  other¬ 
wise  digits  are  dropped  from  tha  right  hand  side  of  the  numbers.  Since  the 
pressures  can  be  expected  to  exceed  one  million  psi,  it  is  advisable  to 
check  the  printed  out-put  against  the  plots. 
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star tc  ru02 

— . — QlM&miQN  jtxnyun.rrimBiiM tTUULMai >m^^.La)jtiui.Lk.LiLj>jLrt^ ul 
101 J »UI 1001 ) »TAU2( 1001  I .YYYI1001 ) iUUI 1001 ) *  T I Tlk  (  26 )  iTITU J) 

& . Uf  Itlil  giUAt&AM  DL-ttEIM  MPKtfWMILMfiLU&kD ..OH..A  WMiiUIfclLJEAillJJI _ 

C  OTHER  THAN  THAT  AT  N.S.R.D.C*  ELIMINATE  THE  NEXT  lAKU 

.T| 


^ii:MJTTi»!?MN*Tjr]WTOKrT^rrMS^<iyTn?PJ|  ’  " 

\fTTTWrW!WifrumtiWTTi 
^■^ltT7nfTjc3Trj«!^T^liuijlTPfn?ii 

yssii  -n*js  i 

vnnlutwTlHufnil 


C  INITIAL  VALUES 

**LlfcN- 1 • 0 ) / 1 2 • 0*tN !  ) 


i‘Li*n 


HO«CN* I PL+B 1 / I OL* (EN-1.0) )*<  ( ( PO+B ) ✓ (PL+flT) *» (< kN- 1 ,0 TV 
*G*PO«MUQ/RO) 


1 OHO /CO 

0'C0",CL*<eN-l.OI*DPO/(2.0*EN*(PL+tt>>*(  (°L+b)/(PO+d)  )**(  (EN+1.0J/2.Q 


1*ENI 


f>DHO»OOPO/OL*  1  <PL+8)/(P0+B)  )»*ll.O/EN)-DPO**2/(EN*DL»(PL+B>  )*(  (PL+ 
) / ( PO+B ) )  *#  (  ( EN+1 ■ 0 ) / kN ) 


.  U*CO  I  /  (CQ-UU  )  +DUNO*  I  DUKO-DCO )  /  (  CO-UO )  ♦ 
UO**2/ 1 2. 0*RO ) * ( DDRO-3 .0*DCO ) / ( CO-UO > +DHO/NO* ( CQ+UO ) / ( CO-UO » +HO/RU 


2*  ( DCO+OORO )  / 1  CO-UO »  +OHO/CO*  ( OCO-ODKO )  /  ( CO-l'O  >  ♦  ( (  UO*DHU+RO#DDMO  )*CU 
/ ( RO*CO**2 ) 


OOCO—CL*  ( EN- 1 . 0  )  »  (  EN+1 .  U I  H)PO»*2/  (  2 . 0*tN*  <  »L+d )  )  **2*  I  (PL+B )  /  I  PO+B 
)  )»»(,(  3. 0*6N+ 1.0)  /(  2.0*kN  >  )+CL«(tN- 1.0 ) »DDPu/ ( 2.0*tN» ( PL+B )  )*<  ( PL+ 


+  1.0)/ (2.0»EN)  I 


WHKmil 

rnrwnsn 


11 .0  J  *UO**3/RO**2-DDORO ) 


wjyrvafr  tvu'u  n'lia'iiu  m  n'Li^iyTiJ  uyn'] 


1) )*( (PL+B ) / ( PO+b ) )  **  ( ( EN  +  1 «  0 ) /EN I +UPO**3* (EN  +  l«0)/( EN#*2*DL* ( PL+d ) 
?**?!*< < PL+B ) / < PO+H 1 »#*( (2.0*EN+1 .0>/fcN) 


OOOORO" ( 6. 0*UO»DDKO**2-4 .0*ODKO**2*CO-3 . 0*UO*OCDRO*CO-« .  0*UO*Ul)Ku* 
/ ( KQ«  <  CO-UO ) )  +  ( 3 .0»UDKO»DUOKU-UURU*DDCQ-2 . Q#DDORO*DCU ) / ( CO-UO ) 


2+13. 5*UO**2*DDDRO- 1 . 5*UO**2*DOCu+ODHU*CO+2 . 0*DHQ#DCO+MC*DDCO ) / <  RO* 
_3  C  CO-UO  ))♦< MO»DDDRO+3.0*DDNO*DHO+3 .0*UO#DDHO+RO*Dl)DHO )  /  (  NO#  (  CO-UO  >  ) 
A 


62 ,0*RO#UO#ODHO»OCO+RO*UO*DHO»DDCO ) / ( KO*CO#*2» ( CO-UO )) -2 . 0*UO*DHO*D 
7CO**2/ !CO##3#  (  CO-UO ) ) 


DD0CO«CL*(EN-l.Q)*( EN+1.0)*<  3.0*fcN+1.0)*DP0##3/(2.Q*EN*(PL+B) )##3* 
II ( PL+B ) / ( PO+B ) )**( ( 5.0*EN+i .01/(2. 0»EN) )-3.0#CL*(EN-1.0)#( EN+1 .0 )# 


0/(2.Q*EN*(PL+B) )**2*( ( PL+B )/(PO+B> )*#(T3iO*EN+ 1.0)/ (2.0* LN 
3) »+CL*(EN-1.0)*DO0PO/(2.O*EtT*(PL+b) )*( ( PL+B )/( PO+b > ) ** ( ( tN+1 .0 ) / ( 2 


* ( 3. 0*0+2. 0 ) #DOKO*UO»»2/NO**2+ 


13.0* ( 3 .0*0+1 .0 ) *OORO»*2/RO-OODONO+3 . 0* ( 3 .0*0+1 .0 ) * ( 3. 0*0+2. 0 ) *  < 0+1 
2.0) *UQ**4/R0»* 3+4.0* ( 3 .0*0+1 .0 ) *UO*ODURO/NO ) 


ODOOHO«ODDOPO/DL*( ( PL+B) /(PO+B) ) ** ( 1 . 0/tN ) - ( 4. 0*UD0P0*DPO+3.0*DDPo 
1**2  >/(EN*DL*< PL+B) )*( ( PL+U ( / < PO+B > )**( ( tN+1 . 0 ) /EN ) +6. 0* ( EN+ 1 . 0 ) *OP 
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20«»2*0L>P0/  (EN»*2»0U*«HL+Bt**2  >«<  i  PL+B  J  /  t PO.d I  )**(  1 2  «0*EN+1.0)  /EN  I* 
..3QPQ****  t  EN+l.Q  I *<2«VI>EN+1.QI/(EN**3*QL<>1PL»B)»43)*I  (PL+ttl / tPo+Bl  )♦ 
Ml  (S.0*eN+1.0)/EW> 

ODPPUO-  I »  I  i.O*PDRQ*PPPRO*CO-12  «Q»UUKu»»2<>P(.0-t «0»Um>DUPURu*Cu-l?  .0 
l»UO»PPPRO*OCtl-12.0»UO*PORO»PDCO-3.0*KO*DPPORO*PCU-J.O*RP«UPPKO*PUC 
20+R0*PPRU»PPPCO4'S  .0*RO«POPHU*»  2+4 .0*KU«UPKO«uPUPHO- 1 .  !>«U0**2*UUUtu 

-^§ifeA368£gijg^8Sg£S^;£giiSK^og£?^3^ggBgS£gOg!i:S:ii:Sii5g:!;>?g8S5 

6*DCQ/CQ»«2+RQ»DPDPHO-9»Q«UU»DPRO»PPrtWCU-3.Q*UP«»2»QPPHQ/CQ*».Q»UP 
7*#2#OOMO*PCO/CO**2-S.O«KO*PPORO*DUMO/CO-J.O*RO*OPKO*pPPHO/COt6.0*R 
-IQ*D3in*DDH0*DCQ/CQ«*2»RPl»ll0*DPDDHU/C-U'»3.Q*RQ*U04PDDHQ*DCl)/CQ-**2+3* 
90*UO**2*PMO*POCO/CO*«2-6.0*UO*#2*PHO*PCO**2/CP**3+3.0*RO*POORO*PMO 
l*DCO/CO4*2+3»0*RL#DPRO*OH0*DPCO/CU»*2'-6  »0*WO*DDRO*DHO*DCO4*2/C04*3 
2'f3.0*RO*UO*PUHO*PPCU/CO**2-6.0*Ru*UO«PPMO«t)CO»*2/CO*»34-HO»UO*PHO*U 
1DQCQ/Cfl**2.-A*Q4RQ*UU*PHP*PIU*DDCU/CQ* «  3 +6 . Q»RP*UO«DHQ«DC0*»3/C0*»» 
4-RO*DDDDRO*DHO /CO ) / <  RO* ( CO-UO )) 

-.ULflfcfl _ _ 

DO  10  \-1.4 

YY(  1  .ll-R0+UO»T+DPRO*T*42/2.0+T«PUPKu»T«*2/6.Q+T4»2»DDPDt<O»T4»2/24 
1.0+T*43*DDDDUO*T««2/120.0 

YY(1 .2)-UO+PORO»T+DDDRO»T»»2/2.0+T*DDDDR04T»*2/6.Q+T»42»PODPUO«T»* 
12/24.0 

YYtl.4|-pp*(RQ/YYUtl))»«(3.Q*Q) 

call  Accel (  Yt (  i  . i  i » yy i 1 . 2 i . y y  <  J . : 


.  2). Vv I  I « 3  >1 YY ( 1.4)) 


ID- 


T=rrf 

Rl-YYIl.l) 

femrnr 


R3-YYI3.il 

1U-VV<4.1> 

Ul-YYI 1 .2  I 

ua-yynr^T 

-YYI3.2I 


# 


U4-YYI4.2 I 
0U1-YY (1.3) 

■w7t;m 


aU2‘ 

DU3-YY 13.3) 
PU4-YYI4.3I 


UP4- 


-U4 


TFT 

C4-U4 


TIT 

0 


04-R4 

1-4 

T^TJT 


0  U*3-Ul ♦4.0*$* ( 2. 0*DU4-DU 3  +2 • 0»PU2 1/3.0 


230 


RP3.Rlj4.Q4S; 1 2 . 0»U4-U3+2 . 0 
COO-OP5-1 12.0*1 UP4-C4 1/121. 


0»U2 1/3.0 


240 


000-RP5- 1 l 2 . 0* I RP4-P4 I / 1 2 1 . 0 

-POT--TOV7T57DC0  H*J*.0*C) - 

CALL  ACCEL ( 000 » COD » OCOP .POD  I 
t5gTT.TnU4.-U243.C«S*rPCOD^2.0»PU4-m3 ) )  /g.O 
03-t9.0*R4-R243.C*S*<CO0+2.04U4-U3»)/«.0 

i  y  rw  i  dps-ci  rrcTTiroTBo 

UPS-C3+9.0* I UP3-C5 1/121.0 

nwswres'g*  <*ps-ds  i  /  nr .  o 


TU  210 


L*1 


C  HALF  INTERVAL 


Y6  230 
PROCEDURE 
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1U2I/120.O 

R3MH«(12»0*R4+135.0*R3+10tt.U**2+Kl >/256«0+S*I-3.0*U4-!>4.0*U3+27.0« 


1U2 )/ 128*0 

U4MH-(80»0*c4n35.0»U3+40.0*U2«-Ul>/256.Q+i>*<-lS.Q*UU4+90.0*DU3  +  15. 


10*DU2»/128.U 

U3HHa( 12.0*U4+133.0*U3+10B.O*U2+Ul)/236.0+S*<-3.0*OU4-34.0*OU3+27. 


10*Dy2l/128.0 


L-L+l 

.1>T.10)WRITEI6»75)  _  _ _ 


75  FORMAT (S8H  THE  PROCESS  IS  NOT  CONVERGING  QUICKLY  ENOUGH  AT  SOME  ST 

FPI 


60  TO  700 

C  CALCULATION  OF  FINAL  VAl  UES 


210  YY(! »2>«C5+9.0*(UPS-C5I/121.0 
00  YY  I ! *  1 ) ■D5+9»0* ( RP5-D5 ) / 12 1 *0 


YY( I «4)-P0*(R0/YY< I ♦ 1 )) ** < 3 .0*G ) 

YYfI»l)tYY(I*2l*YY(I»31»YY(l*4)) 


-1)+S 


C  RELOCATE  CERTAIN  QUANTITIES  FOR  THE  NEXT  STEP  OF  THE  INTEGRATION 

4J! 


04-05 

UP4-UP5 


I'tU'i 


U4-TYI I .2) 
0U1-0U2 


OU2-OU3 


I'lVUH'i 


0U4- YY  (  1 13  i 


7  no  9  J-1»I 

«YY(  J»1I*YY«  J»2l**2/2»0+YY(  J*1  >*{  YY(  J*4)-PL)  /0L«U#0-(  YY(  J*4)-PL) 


1/(2«0»UL«CL**2) I 

•  »2/Y**2)«CL»YV(J*1>«#2*YYU*2MI  L#0-YY(J»2)**2/(2.0»CL*«2I 


l-CL**2*YY< Jtl »/Y*( 1.0-YYI J»2)/CL> 


|iiPiuriMU-Tiiu.i.T4n%niur4 


riiiKi/Lniji 


l»**2*(Y/CL)*»4/I2.0»STF*«4»CL«*2n 

ICL*STF2**2 J* 


1  ♦  <XK/CU**2*< Y/CL>»*4/(2.«STF2**4*CL«*2M 
U-UI J> 


YY  (  J»5)»0L*!Y/STF-U**2/2.0H-DL/t2.Q*CL**2)*<  Y/STF-U**2/2.0)**2 
UaU'J  f  J  ) 


YYY! J)«DL*<  Y/5TF2-U**2/2. )+DL/<2.*CL**2)*nr/STP2-U**2/2. )**2 
TAU!J)"XX! J)+!STF-YV! J.l) >/CL»( I .0-Y Y 1 3  .2 ) »Y Y ! J .1 )/!CL*STF) ) 


TAU2!J)»XX< J)+!STF2-YY< J.l) )/CL*!l.-YY!J.2)*YY(J»l)/!CL*STF2) I 
1 


IF(TAU!J).LT.TAU< J-l) )  1II«J 
.  >AU3tJ).JAU2UJ*l. 

9  XX! J)«XX! J)*1.0E3 

R1TEI6.16 ) H.RO  »P0 .0. B .£N.DLKfc.F»STF.STF2 


16  FORMAT ! 1H1/////////10X.42HGILM0RES  SECOND  ORDER  APPROXIMATION  FOR 
1A  .F6.0.12H  FOOT  WATER  // 1 IX .21HOEPTH  IMPLOSION  OF  A  .F5.1.32H  IN 


2CH  RADIUS  .  *  >.lcU  TO  A  / /?X . 12HPRESSURE  OF  »F<>.2,37H  PSI  WIT 

3H  A  GAS  WHOSE  GAMMA  VALUE  IS  .F6.3//9X.36HB.  N»  AND  DENSITY  OF  THE 


4  LIQUID  ARE  .F8.Q.5H  PSl  .F6.3.5H  AND  //9X.E10.4.32H  LBS-SEC**2/ IN 
5**4  RESPECTIVELY.  .14HSTAN00FFS  ARE  .F6.2.BH  IN  AND  *F6.2.3H  IN) 


IF(III.GT.O)  WRITE! 6. IB) TA  U(Ill) 

18  FORMAT! 1H  //////// 10X . 72H  NON-UN lUUE-VALUES  APPEAR  FOR  STANDOFF  TI 


1MES  IN  REGION  PRECEID ING  TSTF1 ■ . F8. 3 . 10H  MlLLlSECS  ) 

WRITE! 6*11 ) (XX! JI.VY! J.l  ).YY< J.2),YY< J.3) .YY< J»4).T*U!J) .YY! J.5) » 


1U  <  J  ).»  T AU2 !  J  )  *  Y  YY  <  J  )  .UU< J ) . J«1 » I ) 

11  FORMAT! 123H1  TIME  RADIUS  WALL  VEL  ACCEL _ WALL  PRES 

1  T  STFl  P  STF 1  U  STF1  T  STF2  F STF2  CTsT 

2C2  /35H  VALUES  IN  MILLI5ECS*  IPS.  AND  PSI* 


/////!lX.F9*6.F8.5.E13.5.El3.5.fcl2.».FlI.6.El2.4*E12.4. 
3  F10.6.E12.4.E12.4/) ) 


I U fTl  n  4  rTHB;)  J -i- r«j  j  ,-wv ;  I  ■T_W-f«Tr  1  T : 


C  OTHER 

THAN  THAT  AT  N. 

S.R.D.C.  ELIMINATE  ALL  CARDS  BETWEEN  THIS  CmRD  I 

C  AND  THE  NEXT  COMMENT  CARO  WHICH  READS  -  END  OF  PLOTTING  ROUTINE  - 
_ DATA  SCAL.E/6NL1NEAR/ _ 

XL' 
.  CA| 

R.DX) 

DX' 

»2  »0*DX 

i r.v.wi  a  n  frii 


Y2MIN*YY(1.2) 
Y4MAX« YY ( 1 .4) 


Y5MAX«YY< 1.5) 


kfiMiiktuivn-ji 


UMAV.«U<  1) 
UMI*I«U!  1) 


YY ! N . 3 ) »U  <  N ) 


IF!YY!N,2).LT.Y2MAX)GO  TO  8 
Y2MAX»YY ( N»2) 


81  IF(YY(N,2 ) »GT • Y2M IN )GO  TO 
Y2MI N« YY ! N .2 ) 


WVM  Ifffn.fHHl  WCT.TII^B  fill 

IT»  1.11*11:1  ft*  rioj.ii 
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87  1F(UU(N) .LT.UMAX)  GOTO  88 

i 


86  I F ( UU( N  )  .GT »UM IN )  GOTO  80 
(N) 


80  CONTINUE 


|»l»l 

f9a«CT»| 

•  •  '*>£  ■.  j i- * p _ ■ 

iirtTHIrajifl 


READ(5*76)BTITLE( 1 ) »BT I TLE ( 2 ) .bT I TLE I  3 > . bT I TLt ( 4 ) .bT I TLt ( 5 ) *b T I TLt 
) »BT I TLE( 8 ) 


76  FORMAT (10A6) 


DAT  A! T I TLE ( .U  » J- 1 .25 ) /6HRAD 1US *6HI N  IN  *3*6H  .6HVELOCI. 

_ N  .6HINCHES.6H  PER  S»6ntCONL>  »6hWALL  P  .feHKbSSUR  »6HE  IN  Pi 

26HSI  »6H  »6HEULER1 .6HAN  VtL»6hOCITY  .6H1N  IPS»6H  . 

_ 36HQVERPR.6HESSURE.6H  IN  PS»6HI _ jl611 _ L _ 

D071  MI-1.5 

(Mi l*TI TLF (MI  I 


CALL  GPLOTI SCALE. 1. I»0*XL»X**0.©«YT »DX  »  DY ) 
MAX.YT.DYT ) 


Y2MI N--Y2MI N 

CALL  SPACE! Y2MIN.YB.DYB) 


YB— YB 

DY-AMAX1 ( DYB.DYT 


DO 7 2  MI-1.5 

II-MU5 


(MU-TITLE  (11) 

CALL  GPLOT(SCALE.2. I  .0  .XL .XN. Yb. YT .DX .DY ) 


ALL  SPACE(YAMAX.YT.DY) 
D073  MI-1.5 


■MU 

73  TI TLEV (Ml ) -TI TLE ( 1 1  ) 


UMIN— UMIN 

CALL  SPACE (UMIN.YB.DYB) 


DO  55  N-l.I 
YY(N.3)-U( 


55  XX ( N ) -T AU( N ) 


.UMJ.T.VVVr.i'IVVRV.^rlMII 


CALL  SPACE ( A.XR .OX ) 


D074  MI-1.5 


74  TITLEV(MU-TITLE(  II) 


D0135  N-l » I 
YYIN.3 )-UU(N) 


155  XX  ( M ) -TAU2 ( N ) 

CALL  GPLOT(SCALE»3. I.2.XL.XR.Y8. YT .DX .DY ) 


( YSMAx . YT  *DYT ) 


CALL  SPACE! YSMIN.YB.DYB) 


iTTfTTTuni 


196  XX ( N ) -T ACM N ) 
0077  MI-1.9 


0Y-AMAX1 (DYB.DYT ) 
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I I*MI+20 


77  TlTLEV<Mn»TirL£LH J _ 

CALL  GPLOTISCAi-t.S.  i  .1  »XL#XK»YB.  YT.DX.DY  ) 


_ D056N*  1 » I _ _ _ 

YY ( N  #5  )  *  YYY ( N ) 

56  XX I N  )  *  T AU2  ( N  ) _  _ _ 

CALL  GPlOTI SCALE .5. I .2.XL *XK. YB. YT.DX  ,0Y  ) 

C _ END  OF  PLOTTING  ROUTINE _ 

00  TO  22 

999  STOP _ _ _ 

CND 

sibfk  ACCELS _ _ _ 

SUBROUTINE  ACCEL ( R .U .DU . P ) 

_ COMMON/AC/CL »PL »DL .B.EN.G.PU.RO _ 

C*CL«< ( P+8 ) / <  PL+B ) )**< IEN~1.0)/<2.Q*EN)  ) 

_ H*E.“1*1  PL+fl  J  /  I  l  EN— 1 »  U )  *DL  I  *  (  (  (  P  +  U  )  /  1  PL+B  )  I  **  ((  EN- 1 .0) /tN ) -1 . 0 ) 

DP»-3.0#G*P0*U/R*<R0/R)#rf(3.0*G) 

DH«DP/PL»t lPL+B)/tP*B)  )»»n  .O/ENl _ 

DU*U**2/ ( 2.0*R  >*< U-3.0*C ) / ( C-U ) +H* ( C+U ) / ( R»  <  C-U i ) +UH/C 


RETURN 

END 


J I  BF.T.C  .-SPACES 


NODEOC.SDD 


SUBROUTINE  SPACE ( ENDPT . XR »DX } 


IF ( UNDPT. Lh« 1 «QE-3 )G0  TO  66 
X  I  D-*ALOG10<  ENDPT*  1000.0) -3.0 


1 F  (  X  I D  .  GE  .  0  -  0  / GO  TO  114 
ID*X ID-1.0 


Atf  .Tfi  61 


114  ID-XID 

1 F ( I D.GE.l ) GO  TO  366 
61  I XR* ENDPT* 10.0** ( 1-ID) 


_ XR*UR±i _ 

XR-XR/10.0**! 1-ID) 

_ If  1IXR..LE,20>GO  T9  Ml 

DX*.' XR/20 
DX*DX/1Q.0»*( 1 —  I D I 
GO  TO  111 

361  DX»1.0/lO.O**( 1-ID) 

GO  TO  111 

366  IXR*ENPPT/10.0*«( ID-1  ) 
XR* I XR+ 1 


AaraR»i£ifrvu-p-i) _ 

IF( IXR.LE.20) GO  TO  136 

DX* I XR/20 _ 

DX»DX*10.0**( ID-1 > 


J&fl  TO  in 


136  DX* 1 .0* 10.0** ( ID- 

3Q  IQ  at 


1  ) 


68  WRITE  I  6. 67) 

67  cORMAT( 1H1///1QX.46H 


TMfc  tNUPOlNT 


1  ) 

111  RETURN 


SIBLDR 

mur 


END 

MOLLV 

TCW 


S1BLDR  VCHARV 
SIBLDR  TABL1Q 
SIBLDR  RITE2Q 
S I  BLDtt  APflNTV 
SIBLDR  APLOTV 


IS  TOO  SMALL  F UW  THIS  SUBKU uTlNt 


HOLLOOOO 

PloToCoc" 


_ __  VCHAOOOQ 

10/10/65  T  AULOUUU 

_ _____  HITtOUOU 

6/15/65  APHNOOOO 
10/04/65  AFLUOOoO 


0*/l0/*fc 


wPiuuuuu 


SlBlDH  OPLOT 
_ JU241A  . . 


Hi 


For  each  bubble  collapse,  Information  must  be  read  in  on  data  cards 
in  tho  following  wuy: 


Card 

1 

Cols 

1-12 

collapse  depth  in  feet  of  water, 

F12.S 

Cols 

13-24 

initial  sphere  radius  in  inches, 

F12.5 

Cols 

25-36 

initial  internal  gas  pressure  in  psi 

F12.S 

Cols 

37-48 

the  specific  heat  ratio  Y  for  the  gas 
inside  the  sphere. 

F12.& 

Cols 

49-60 

#1  standoff  in  inches, 

P12.S 

Cols 

61-72 

*2  standoff  in  inches 

F12.5 

Cols 

73-80 

the  value  of  n,  the  exponent  in  the  equation  of 
state  for  the  adiabatic  compression  of  the 
liquid 

F  8.5 

Card 

2 

Cols 

1-15 

the  value  of  B  in  psi,  a  constant  in  the  equation 
of  state  for  the  adiabatic  compression  of  the 
liquid, 

E15.5 

Cols 

16-30 

the  density  of  the  liquid  in  lb-sec2/in*  at 
standard  temperature  and  pressure, 

E15.5 

Cols 

31-45 

the  sound  speed  in  the  liquid  in  in/sec  at 
standard  temperature  and  pressure. 

E15.5 

Card 

3 

contains  graph  labels.  Jf  the  program  is  not  being 

used 

at  NSRDC,  then  this  card  must  not  be  included  in  the  data. 

Card  3 

Cols  1-48  Main  graph  title  for  all  graphs. 

Two  blank  cards  in  succession  stop  the  computer.  Each  case  requires  no 
more  than  3  minutes  running  ime. 
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Two  methods  are  presented  for  calculating  the  instantaneous 
pressure,  velocity,  acceleration,  and  radius  associated  with  the 
collapse  of  a  spherical  gas-filled  cavity  in  an  infinite  com¬ 
pressible  liquid.  One  is  based  on  the  ideal  gas  law,  the  other  is 
based  on  the  Beattie-Bridgeman  equation  of  state  for  the  gas  inside  the 
cavity.  In  most  cases  the  latter  assumption  must  be  restricted  to 
relatively  mild  implosions.  The  good  agreement  between  the  two  methods 
serves  to  verify  their  validity. 

Included  are  listings  of  the  two  Fortran  IV  computer  programs 
used  to  obtain  numerical  results  of  the  analyses  based  on  the  ideal 
and  Beattie-Bridgeman  gas  models.  The  influence  of  several 
different  gases,  initial  internal  gas  pressures,  and  liquids  on  the 
collapse  is  studied.  On  the  basis  of  explanations  of  the  resulting 
behavior,  new  methods  of  producing  similar  behavior  are  discussed. 
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